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Summary 

r' 

A  code  for  computing  the  aerodynamic  design  of  a 
multistage  axial-flow  compressor  and,  if  desired,  the 
associated  blading  geometry  input  for  internal  flow 
analysis  codes  is  presented.  The  aerodynamic 
solution  gives  velocity  diagrams  on  selected 
streamlines  of  revolution  at  the  blade  row  edges. 
Blading  is  defined  from  stacked  blade  elements 
associated  with  the  selected  streamlines.  The  blade 
element  inlet  and  outlet  angles  are  established 
through  empirical  incidence  and  deviation  angle 
adjustments  to  the  relative  flow  angles  of  the  velocity 
diagrams.  The  blade  element  centerline  is  composed 
of  two  segments  tangentially  joined  at  a  transition 
point.  The  local  blade  angle  variation  of  each 
segment  can  be  specified  with  a  fourth-degree 
polynomial  function  of  path  distance.  Blade  element 
thickness  also  can  be  specified  with  fourth-degree 
polynomial  functions  of  path  distance  from  the 
maximum  thickness  point. 

Steady  axisymmetric  tew  js  assumed;  so  the 
aerodynamic  problem  car^fi^  reduced  to  solving  the 
two-dimensional  flow  field  in  the  meridional  plane. 
Because  the  equations  of  motion  as  developed  herein 
are  only  applicable  for  calculation  stations  outside 
the  blade  rows,  stations  at  the  blade  edges,  but  not 
inside  the  blade  rows,  are  used.  The  streamline 
curvature  method  is  used  for  the  iterative 
aerodynamic  solution.  If  a  blade  design  is  desired, 
the  blade  elements  are  defined  and  stacked  within  the 
aerodynamic  solution  iteration.  Thus  the  design 
velocity  diagrams  can  be  located  at  the  blade  edges. 

The  program  input  includes  the  annulus  profile, 
the  overall  compiessor  mass  flow,  the  pressure  ratio, 
and  the  rotative  speed.  A  number  of  parameters  are 
input  to  specify  and  control  the  blade  row 
aerodynamics  and  geometry.  There  are  numerous 
options  for  controlling  the  way  information  is  input 
and  tor  specifying  the  amount  of  output.  The  output 
from  the  aerodynamic  solution  has  an  overall  blade 
row  and  compressor  performance  summary  followed 
by  blade  element  parameters  for  the  individual  blade 
rows.  It  desired,  blade  coordinates  in  the  streamwise 
direction  toi  internal  How  analysis  codes  and/or 
coordinates  on  plane  sections  through  blades  for 
lubrication  Jiawings  can  be  printed  and  punched. 


Introduction 

Ihc  axial-flow  compressor  is  used  for  aircraft 
engine's  bemuse  li  has  distinct  configuration  and 


performance  advantages  over  other  compressor 
types,  but  the  good  poten'ial  performance  is  not 
easily  attained.  The  problem  and  challenge  to  the 
designer  is  to  model  the  actual  flows  well  enough  to 
adequately  predict  aerodynamic  performance. 
Progress  is  continually  being  made  with  codes  for 
computing  the  complex  three-dimensional  flows  in 
turbomachinery.  However,  it  is  extremely  difficult  to 
design  mechanically  acceptable  turbomachinery 
blading  by  using  the  direct  approach  (i.e.,  specifying 
inviscid  blade  surface  velocities  and  computing  the 
blade  geometry).  Consequently,  the  more  detailed 
codes  are  generally  used  in  the  analysis  mode;  that  is, 
the  flow  field  is  calculated  for  a  fixed  geometry.  The 
current  procedure  is  to  establish  blading  geometry 
with  simpler  design  codes  and  then  to  use  the  more 
detailed  analysis  codes  in  blade  rows  where 
troublesome  flow  conditions  are  likely  to  exist.  In 
this  way  prototype  designs  can  often  be  adjusted 
before  hardware  is  built  and  tested. 

The  time  and  effort  needed  to  get  acceptable 
configurations  can  be  reduced  if  the  design  code  can 
be  made  to  yield  a  good  initial  solution  and  if  the 
design  and  analysis  codes  can  be  made  more 
compatible  with  one  another.  This  compatibility  can 
be  achieved  (1)  if  the  output  from  a  design  code  can 
be  directly  used  by  flow  and  mechanical  analysis 
codes  and  (2)  if  corrective  adjustments  indicated  by 
the  analysis  codes  can  effectively  be  made  in  the 
design  code.  With  these  objectives  in  mind  a 
composite  aerodynamic  and  blade  design  code  for 
axial-flow  compressors  has  been  developed.  The  code 
and  its  capabilities  are  the  subjects  of  this  report. 

The  aerodynamic  solution  assumes  steady, 
axisymmetric  How  and  uses  a  streamline  curvature 
method  for  calculation  stations  outside  the  blade 
rows.  The  program  is  structured  so  that  the  empirical 
correlations  (such  as  those  for  loss,  deviation  angle, 
and  incidence  angle)  can  readily  he  changed  when  the 
need  or  desire  exists.  The  method  of  describing 
blading  is  a  compromise  between  the  vast  amount  of 
input  needed  for  completely  general  blade  elements 
and  the  restrictions  of  simple  shapes.  A  blade 
element  is  defined  on  a  conic  surface  with  thickness 
applied  us  a  centerline  that  is  composed  of  two 
segments  tangentially  joined  at  a  transition  point 
The  blade  angle  function  of  each  segment  can  be 
defined  with  a  fourth-degree  polynomial.  Thickness 
is  prescribed  by  f  irst  specifying  a  maximum  thickness 
value  and  location.  The  distribution  of  thickness  in 
each  direction  from  the  maximum  thickness  location 
is  then  prescribed  with  a  fourth-degree  polynomial. 
Finally  each  polynomial  coefficient  is  defined  across 


blade  elements  with  a  third-degree  polynomial 
function  of  annulus  height. 


Compressor  Design  Procedures 

The  discussion  of  the  compressor  design 
procedures  is  organized  according  to  usage  in  the 
computer  program;  so  for  better  orientation  an 
operational  overview  of  the  program  is  given  first 
(table  1).  The  computer  program  can  be  divided  into 
three  major  phases  of  calculation:  (1)  the  input  and 
initialization  phase,  (2)  the  iteration  phase,  and  (3) 
the  terminal  calculation  phase.  In  the  input  and 
initialization  phase  the  input  data  are  read  and 
interpreted,  the  calculation  stations  are  located  with 
estimated  values  for  the  blade  edges,  and  streamlines 
are  located  on  the  basis  of  annulus  area.  Estimates  of 
stagnation  temperature  and  pressure  and  axial  and 
tangential  velocity  components  are  also  made  for  all 
calculation  points  in  the  flow  field. 

The  iteration  phase  includes  both  the  flow  field 
and  the  blade  design  iterations.  In  the  flow  field 
iteration  the  equations  of  motion  are  satisfied  in  the 
meridional  ( r-z )  plane  for  stations  that  are  lines 
across  the  flow  annulus.  At  the  stations  the  equations 
of  motion  and  overall  flow  continuity  are  satisfied 
with  fixed  values  of  streamline  slope  and  curvature 
for  a  complete  computational  pass  across  the 
annulus.  After  the  overall  flow  continuity  condition 
at  a  calculation  station  is  satisfied,  the  internal 
streamline  intersections  with  the  station  lines  are 
updated  by  solving  for  the  locations  that  give 
specified  fractions  of  overall  station  weight  flow.  At 
the  completion  of  a  pass  through  all  the  calculation 
stations  in  the  annulus,  the  new  streamline  locations 
are  curve  fit  for  new  streamline  slope  and  curvature 
values. 

To  insure  proper  location  of  the  blade  edge 
stations,  most  of  the  blade  design  iteration  is  made 
concurrent  with  the  flow  field  iteration.  This 
operation  includes  the  calculation  of  incidence  and 


deviation  angles,  the  layout  and  stacking  of  blade 
elements,  and  the  realignment  of  the  elements. 

The  terminal  calculation  phase  performs  the  final 
calculations  and  generates  the  output.  Mass-averaged 
parameters  for  the  individual  and  cumulative 
compressor  blade  rows  are  computed  and  printed 
first.  Then  tabulated  values  of  aerodynamic  and 
blading  parameters  along  the  station  lines  are 
computed  and  printed.  Finally  blade  section 
coordinates  and  other  section  mechanical  properties 
can  be  computed  and  printed  if  desired. 

The  program  is  discussed  in  greater  detail  in  the 
following  subsections. 

Input  and  Initialization 

The  basic  computational  plane  is  the  meriodional 
(r-z)  plane  of  a  cylindrical  coordinate  system.  A 
graphic  view  of  an  example  compressor 
configuration  is  shown  in  figure  1 .  The  hub  and  tip 
casing  walls  are  fixed  input.  Calculation  stations  are 
located  at  the  blade  row  leading  and  trailing  edges 
and  at  other  annular  locations  for  the  purpose  of 
locating  streamlines.  The  input  data  can  be  classified 
into  two  groups:  general  information  and  calculation 
station  and  blade  row  information.  The  input 
parameters  and  options,  along  with  the  input  data 
format,  are  described  in  appendix  B.  (All 
mathematical  symbols  are  defined  in  appendix  A.) 
Additional  advice  on  how  to  set  up  the  input  is  given 
in  the  section  User  Information. 

General  Information 

All  the  general  information  is  read  in  first. 
Included  are  the  following: 

(1)  Compressor  rotational  speed 

(2)  Inlet  flow  rate 

(3)  Desired  compressor  pressure  ratio 

(4)  Gas  molecular  weight 

(5)  Number  of  streamlines 
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Figure  1.  -  Calculation  stations  in  compressor  flow  path. 


(6)  Number  of  blade  rows 

(7)  Number  of  annular  stations 

(8)  Coefficients  for  cp  as  a  fifth-degree 
polynomial  function  of  temperature 

(9)  Far  upstream  values  of  total  temperature, 
total  pressure,  and  inlet  tangential  velocity  for  each 
streamline 

(10)  Streamtube  mass  flow  fractions  between 
streamlines 

(11)  Sets  of  points  to  define  tip  and  hub  casing 
contours 

(12)  Sets  of  blade  element  profile  loss  parameters 
that  are  tabulated  as  functions  of  blade  element 
loading  parameter  and  fraction  of  passage  height 

As  many  as  five  loss  sets  can  be  input.  The  particular 
loss  set  used  for  a  given  blade  row  is  designated  in  the 
blade  row  input.  Usually  at  least  two  loss  sets  are 
input — one  for  rotors  and  another  for  stators. 

Calculation  Station  Data  Sets 

The  data  sets  that  contain  information  about  the 
calculation  stations  and  blade  rows  are  read  in  order 
from  annulus  inlet  to  outlet.  The  first  card  of  the 
data  set  identifies  the  type  of  station,  the  tip  and  hub 
axial  locations,  the  tip  and  hub  boundary  layer 
blockage  factors  and  the  station  mass  flow  bleed.  For 
annular  stations  the  single  card  is  the  whole  data  set. 
For  rotors  and  stators  several  cards  are  used  to 
describe  (1)  the  blade  row  inlet  and  outlet  station 
information,  (2)  the  blade  row  aerodynamic 
parameter  input  and  controls,  and  (3)  the  parameters 
defining  blade  geometry.  A  blade  and  the  associated 
edge  calculation  stations  are  located  in  the  annulus 
by  using  a  reference  blade  element  stacking  line. 
Stacking  axis  tip  and  hub  axial  locations  and  lean 
angle  in  the  circumferential  direction  are  input. 

The  locations  of  the  calculation  stations  at  the 
blade  edges  are  at  first  approximated  from  some  of 
the  input  blade  geometry  information.  The  station 
locations  are  moved  during  later  iterations  when  the 
blade  elements  are  defined  and  stacked.  However, 
the  input  tip  and  hub  boundary  layer  blockages  and 
mass  flow  bleeds  for  the  inlet  and  outlet  stations  are 
constant. 

Aerodynamic  parameter  input  and  controls. — The 
blade  aerodynamic  design  is  controlled  with  several 
parameters  that  impose  the  necessary  and  sufficient 
conditions  for  a  solution.  The  options  as  to  how  such 
conditions  can  be  imposed  are  shown  in  table  II.  For 
rotors  the  most  convenient  option  is  to  specify  the 
stage  energy  addition  as  a  cumulative  fraction  of  the 
overall  compressor  energy  addition.  With  this  option 
the  radial  distribution  of  energy  addition  is  not  input 
directly  but  is  imposed  through  a  normalized  rotor 
exit  stagnation  pressure  profile  that  is  expressed  as  a 


polynomial  function  of  annulus  height  in  the  radial 
direction.  The  pressure  level  is  computed  internally 
to  the  program  from  the  energy  input  level  and  the 
computed  losses.  With  the  other  rotor  options  the 
exit  temperature  profile  is  input  instead  of  the  energy 
addition  fraction  being  specified.  For  either  a  rotor 
or  stator,  stagnation  pressure  profiles  can  be  input 
instead  of  the  losses  being  computed  internal  to  the 
program.  These  options  can  be  useful  to  users  who 
have  existing  aerodynamic  designs  but  want  to  use 
this  program  for  blade  description  and  fabrication 
coordinates. 

At  a  stator  exit  a  tangential  velocity  profile  is  input 
as  a  polynomial  function  of  radius.  Unless  specified, 
the  stator  outlet  pressure  profile  is  determined  from 
stator  losses  and  streamline  mixing  effects  from  the 
upstream  station. 

There  are  some  input  aerodynamic  limits  that  the 
program  will  not  allow  to  be  exceeded.  For  a  rotor 
the  limiting  parameters  are  tip  diffusion  factor  and 
absolute  flow  angle  at  the  hub.  The  stator 
aerodynamic  limits  are  diffusion  factor  and  inlet 
Mach  number  at  the  hub.  If  an  aerodynamic  limit  is 
exceeded  during  iteration,  the  stage  energy  addition 
is  lowered  by  the  amount  needed  to  get  the 
aerodynamic  limit  within  bounds.  If  any  other  stage 
is  not  up  to  one  of  the  aerodynamic  limits,  the  energy 
decrement  is  made  up  among  such  stages.  If  all  the 
stages  reach  an  aerodynamic  limit,  the  input  overall 
compressor  pressure  ratio  is  lowered. 

The  blade  angles  are  related  to  fluid  flow  angles 
along  streamlines  by  two  key  correction 
parameters — incidence  angle  at  the  inlet  and 
deviation  angle  at  the  outlet  (fig.  2).  There  are  several 
options  for  specifying  each.  Two  of  the  options  for 
both  the  incidence  and  deviation  angles  are  the  two- 
and  three-dimensional  methods  of  reference  1.  Tl^e 
other  incidence  angle  option  is  user-entered  tabular 


Figure  2.  -  Blade  element  incidence  and  deviation  angles. 


data  referenced  to  either  the  centerline  or  the  suction- 
surface  blade  angles  at  the  inlet.  Other  deviation 
options  are  user-entered  tabular  data  and  a  version  of 
Carter’s  rule,  which  was  modified  to  account  for 
centerline  shapes  other  than  a  circular  arc.  The 
modification  is  shown  in  figure  3. 

Another  input  aerodynamic  parameter  is  the 
minimum  blade  choke  margin  (A  I A  *)  - 1 ,  where  A  is 
the  local  streamtube  cascade  channel  area  and  A*  is 
the  corresponding  area  needed  for  choked  flow.  The 
A *  value  is  the  area  needed  to  pass  the  streamtube 
flow  at  a  relative  Mach  number  of  1 .0.  The  effects  of 
losses  in  all  blade  rows  and  energy  addition  in  rotors 
are  included  in  the  computation  of  A*.  Choke 
margin  depends  on  the  flow  conditions  and  geometry 
defining  the  channel  area.  If  insufficient  choke 
margin  exists  in  a  prototype  design,  some 
compromise  must  be  made  in  either  the  aerodynamic 
requirements  or  the  geometry.  Minor  choke  margin 
deficiencies  can  usually  be  accommodated  with 
adjustments  in  geometry.  Logical  procedures  for 
geometry  adjustments  are  not  obvious;  however,  if 
the  minimum  margin  occurs  at  the  channel  entrance, 
increased  incidence  is  an  effective  method  of  relief.  If 
a  minimum  desired  choke  margin  is  input,  the 
program  will  adjust  incidence  angle  up  to  +  2°  to  the 
leading-edge  suction  surface  in  order  to  attain  the 
specified  choke  margin  if  the  channel  entrance  is  the 
problem.  When  the  minimum  margin  occurs  at  other 
locations  in  the  channel,  the  minimum  value  and  its 
location  are  printed  in  the  output  and  it  is  up  to  the 
user  to  decide  if  he  wants  to  make  compromises  to 
improve  the  choke  margin. 


Figure  3.  -  Variation  of  coefficient  for  Carter's  deviation  equation 
with  location  of  Wade  element  maximum  camber  point, 
m  - 10. 219  ♦  0. 0008916  r  *  0. 000027085 
x(2a/ci  (2.  175-0. 035525  r*0. 00010167^1 


Blade  geometry  parameters. — A  number  of  blade 
geometry  parameters  are  input  for  the  pupose  of 
defining  a  blade.  Blade  chord  is  defined  along  flow 
streamlines,  but  for  the  purpose  of  this  blade 
definition  a  radial  projection  of  streamline  chord  is 
specified  because  it  is  more  meaningful  for  defining  a 
structurally  sound  configuration.  The  radially 
projected  chord  is  defined  from  the  number  of 
blades,  the  tip  solidity,  and  a  normalized  polynomial 
for  the  radial  variation  of  chord.  The  blade  is 
basically  defined  from  a  stacked  series  of  gradually 
changing  airfoil  shapes  or  “blade  elements”  in  the 
radial  direction. 

Each  blade  element,  as  shown  in  figure  4,  is 
defined  from  a  thickness  distribution  applied  to  a 
two-segment  centerline.  The  variation  of  the  local 
centerline  angles  k  with  path  distance  can  be  specified 
by  option  through  the  parameter  IDEF(IROW).  If 
IDEF(IROW)  equals  zero,  the  k  for  each  segment 
varies  linearly  with  path  distance  (as  a  circular  arc). 
When  IDEF(IROW)  does  not  equal  zero,  the  k  for 
each  segment  is  expressed  as  a  fourth-degree 
polynomial  function  of  path  distance.  The  blade 
angle  is  continuous  at  the  transition  point,  but  the 
rate  at  which  the  angle  changes  with  distance 
(curvature)  can  be  discontinuous.  The  ratio  of 
curvature  for  the  first  segment  to  that  for  the  second 
segment  is  defined  as  the  turning  rate  ratio.  When  the 
blade  local  centerline  angle  k  is  specified  by 
polynomial  coefficients,  the  turning  rate  ratio  is 
controlled  by  the  relative  magnitudes  of  the  linear 
term  coefficients  of  the  polynomials  for  each 
segment.  However,  when  the  segments  are  treated 


Figure  4.  -  Reference  and  direction  nomenclature  for  prescribed  blade  element 
centerline  and  thickness  polynomials. 


simply  as  circular  arcs,  the  turning  rate  ratio  is  a 
blade  element  input  parameter. 

When  lDEF(IROW)  equals  zero,  there  are  some 
options  for  specifying  the  turning  rate  ratio  at  the 
transition  point.  With  the  CIRCULAR  option  the 
value  is  set  at  1 .0,  as  for  a  circular  arc  blade  element, 
for  all  blade  elements  in  the  blade  row.  With  the 
TABULAR  option  a  table  of  values  for  the  elements 
is  read.  With  the  OPTIMUM  option  a  value  will  be 
set  by  an  empirical  function  of  inlet  relative  Mach 
number.  For  this  option  the  blade  element  will  be  a 
circular  arc  below  a  relative  Mach  number  of  0.8.  As 
relative  Mach  number  increases,  the  ratio  of  first-  to 
second-segment  turning  rate  at  the  transition  point  is 
reduced.  A  limit  of  zero  camber  on  the  suction 
surface  of  the  first  segment  is  approached  at  an  inlet 
relative  Mach  number  of  about  1.60. 

The  coefficients  for  the  centerline  polynomial  (i.e. , 
when  IDEF(IROW)  ^0)  are  input  as  a  cubic  function 
of  blade  span.  There  are  two  reasons  for  this  method 
of  specification.  First,  the  user  is  more  confident  of 
specifying  a  relatively  smooth  blade  surface;  and 
second,  the  amount  of  input  is  reduced  over  that 
required  by  individual  coefficients  for  as  many  as  1 1 
blade  elements. 

Blade  element  surface  definition  begins  with  three 
anchor  points  from  the  centerline.  These  points  are  a 
maximum  thickness  point  and  the  two  end  points.  A 
maximum  thickness  value  normalized  by  chord  and 
its  location  as  a  fraction  of  chord  are  input.  At  the 
maximum  thickness  point  the  normal-to-centerline 
distance  to  each  surface  is  one-half  the  maximum 
thickness,  and  the  surface  k  angles  are  equal  to  the 
centerline  *c. 

At  the  blade  element  ends  the  leading-  and  trailing- 
edge  end  circle  radii  normalized  to  chord  are  input.  If 
IDEF(IROW)  does  not  equal  zero,  the  end 
configurations  are  ellipses  with  semimajor  axes 
tangent  to  the  local  centerline.  For  this  case  the  input 
end  circle  radius  is  used  as  the  minimum  radius  value 
of  the  ellipse.  For  each  ellipse  one  other  parameter  is 
input  to  specify  elongation.  The  parameter  is 
e=(b/a)~  1,  where  b  and  a  are  the  semimajor  and 
semiminor  axes,  respectively.  Note  that  as  e 
approaches  zero,  the  ellipse  approaches  a  circle  with 
the  input  radius. 

A  surface  definition  criterion  is  that  the  surface 
curve  join  the  end  circles  or  ellipses  at  a  point  of 
tangency.  When  IDEF(IROW)  equals  zero,  the 
surface  curves  are  defined  with  x  being  a  linear 
function  of  path  distance  for  each  segment.  As 
explained  in  reference  2,  necessary  and  sufficient 
conditions  exist  to  completely  define  the  surfaces 
when  the  computation  is  begun  on  the  segment  where 
the  maximum  thickness  occurs. 

When  IDEF(IROW)  does  not  equal  zero,  the  blade 


surfaces  for  each  segment  are  defined  by  polynomial 
distributions  of  the  normal-to-centerline  distance. 
The  functional  relation  for  this  distance  is 

t=‘-f  -aVs0+a\IS^S0+  ~  -bS2-cS*-dS* 

2vS0 

where  S  is  the  centerline  distance  (normalized  by 
chord)  from  the  maximum  thickness  point.  Values  of 
S  are  positive  in  either  direction  from  the  maximum 
thickness  point;  and  S0  is  the  maximum  S,  which  is 
the  distance  to  the  point  where  the  end  ellipse 
intersects  the  centerline  (fig.  4). 

There  are  two  other  input  parameters  for  blade 
rows.  One  is  a  material  density  for  rotors.  If  a 
nonzero  value  is  input  for  a  rotor,  the  stacked  blade 
will  lean  in  both  the  meridional  and  the  r-6  planes  so 
that  the  centrifugal  force  on  a  blade  with  the  input 
material  density  will  balance  the  aerodynamic  forces 
at  the  design  point.  The  objective  is  to  minimize  the 
blade  root  stress.  With  atmospheric  air  as  the 
working  fluid,  the  lean  is  normally  only  a  fraction  of 
a  degree. 

The  final  input  parameter,  NXCUT,  controls  the 
number  and  location  of  planes  through  a  blade  row 
for  which  fabrication  coordinates  are  desired.  If  the 
parameter  is  zero,  the  program  will  set  the  number  of 
XCUT’s  on  the  basis  of  aspect  ratio,  which  is  the 
ratio  of  overall  radial  to  axial  blade  lengths.  For 
positive  parameter  input  values  the  program  will 
determine  appropriate  locations  for  that  number  of 
planes  to  represent  the  blade.  Negative  parameter 
values  trigger  an  option  to  read  cards  for  the  XCUT 
plane  values.  The  number  of  input  values  expected 
for  a  blade  row  is  the  absolute  value  of  the  negative 
parameter. 

Initialization 

Once  the  input  is  read,  a  number  of  initialization 
calculations  are  made  in  subroutine  START  in 
preparation  for  the  iterative  phase  of  computation. 
The  axial  locations  of  the  blade  edges  are 
approximated  and  the  intersections  of  all  station  lines 
with  the  casing  walls  are  determined.  Checks  are 
made  to  be  certain  that  the  spacing  of  calculation 
stations  is  appropriate.  Annular  stations  will  be 
shifted  by  the  program  if  calculation  stations  cross 
one  another  or  if  adjacent  spacing  is  less  than  30 
percent  of  the  spacing  of  neighboring  stations. 

Streamlines  are  initially  positioned  by  applying 
the  input  stream-tube  weight  flow  fractions  to  the 
annulus  area.  From  the  input  data  the  cir¬ 
cumferential  component  of  velocity  and  the 
stagnation  temperature  and  pressure  are 
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approximated  for  all  streamlines  at  all  calculation 
stations  throughout  the  flow  field.  Finally  an  axial 
velocity  is  computed  for  each  station  by  using 
meanline  values  in  a  continuity  calculation  at  the 
station. 

Iteration 

The  general  objective  of  the  program  is  to  obtain 
both  an  aerodynamic  solution  and  a  blade  design. 
Both  are  achieved  with  iterative  procedures.  The 
aerodynamic  design  has  the  greater  sensitivity,  and  it 
requires  more  iterations.  The  program  is  set  up  to  do 
the  aerodynamic  and  blade  design  iterations 
concurrently.  However,  the  blade  design  is  done  less 
frequently  and  lags  the  aerodynamic  iteration.  The 
first  blade  design  iteration  occurs  on  the  fourth 
aerodynamic  iteration,  and  the  final  blade  design 
pass  is  made  after  the  aerodynamic  solution  is 
printed. 


Aerodynamic  Design 

The  aerodynamic  design  solution  establishes 
complete  velocity  diagrams  and  fluid  state  properties 
on  streamlines  at  the  blade  row  inlet  and  exit.  A 
bilevel  iteration  is  used  to  arrive  at  the  solution.  In 
the  outer  loop  the  variables  are  stagnation 
temperature  and  pressure;  the  tangential  component 
of  velocity;  and  the  streamline  location,  slope,  and 
curvature.  The  inner  loop  is  the  station  flow 
continuity  calculation  in  which  the  axial  component 
of  velocity  is  the  variable  and  the  outer  loop 
parameters  are  held  fixed.  An  example  flow  field 
with  typical  placement  of  calculation  stations  is 
shown  in  figure  1 . 

Outer  loop. — In  the  program  the  control  routine 
for  the  outer  loop  is  VDIAG.  The  basic  procedure  is 
station  marching  from  inlet  to  outlet  with  streamline 
parameters  fixed.  Only  after  a  pass  through  all  the 
stations  are  the  streamlines  relocated  from  the 
current  flow  solution.  Normally  between  10  and  20 
of  the  cycles  are  needed  to  converge  to  a  solution. 

The  major  part  of  the  blade  design  is  also 
controlled  in  the  outer  loop.  When  a  blade  design 
iteration  is  made,  the  blade  edge  station  locations  are 
moved  to  the  new  blade  edge  locations. 

The  tangential  velocity  and  the  stagnation 
temperature  and  pressure  at  a  station  are  determined 
as  changes  from  values  of  the  preceding  station  on 
the  particular  streamline.  For  annular  stations  and 
blade  row  inlets  the  tangential  velocity  is  determined 
from  the  conservation  of  angular  momentum;  that  is, 
the  product  of  radius  and  tangential  velocity  remains 
the  same  along  streamlines  outside  the  blade  rows. 
Stagnation  temperature  and  pressure  should  also  be 


conserved  along  streamlines  outside  the  blade  rows 
except  for  mixing  effects  from  turbulence  and 
secondary  flows.  The  stagnation  pressure 
distribution  is  input  behind  the  rotors;  so  pressure 
gradients  are  reasonably  well  controlled  in  the  design 
process  without  using  empirical  mixing  terms. 

In  the  design  process  the  rotor  energy  addition 
must  cover  nonproductive  losses  in  addition  to 
producing  a  desired  pressure.  With  the  usual  input 
options,  losses  are  computed  internal  to  the  program. 
Normally  there  is  a  significant  radial  gradient  of  loss; 
so  there  is  also  a  radial  gradient  of  work.  The 
stagnation  temperature  increase  along  a  streamline  is 
in  almost  direct  proportion  to  the  blade  element 
work;  so  temperature  gradients  are  generated. 
Because  these  gradients  through  compressor  stages 
are  basically  additive,  theoretically  the  gradients  can 
grow  very  large.  The  real  flows  in  compressors 
reduce  this  effect  somewhat  with  fluid  mixing.  To  at 
least  partially  account  for  fluid  mixing  in  an 
empirical  manner,  a  mixing  term  for  temperature  is 
used  in  the  program.  The  mass  average  temperature 
is  held  constant  at  a  station,  but  specific  streamline 
values  outside  the  blade  rows  are  modified  from  the 
previous  station  values  by  equation  (1). 


( 3f);-( p  [-"■002(S’)/-iH <» 


where  A z  is  the  axial  distance  between  the  adjacent 
stations.  Future  adjustments  in  this  functional 
relation  are  probable  as  data  from  multistage 
compressors  become  available. 

Stagnation  temperature  and  pressure  values  are  the 
most  difficult  to  set  at  blade  row  exits.  This  is  mainly 
because  of  the  complex  real  flow  effects  through  a 
blade  row  that  must  be  represented  either  through 
theoretical  models  of  loss  or  by  empirical 
correlations.  Representation  of  losses  is,  of  course, 
one  of  the  major  problems  for  an  aerodynamic 
solution.  In  this  program  the  losses  are  represented 
by  two  additive  components:  shock  losses,  and  all 
other  losses. 

The  shock  losses  are  a  modification  of  those  given 
in  reference  3.  This  reference,  in  essence,  gives  the 
shock  loss  associated  with  a  normal  shock  with  an 
approach  Mach  number  equal  to  the  average  relative 
Mach  number  at  the  suction  and  pressure  surfaces  of 
the  blade  at  the  normal  shock.  The  suction-surface 
Mach  number  at  the  shock  is  determined  by  Prandtl- 
Meyer  turning  from  the  inlet. 

Unless  the  flow  is  in  the  low  transonic  range,  a 
normal  shock  cannot  be  maintained  in  a  blade 
channel.  Either  the  shock  is  oblique  or  it  develops  a 
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foot  at  the  blade  surface  because  the  boundary  layer 
cannot  sustain  the  sudden  static  pressure  rise.  In 
either  case  the  shock  losses  are  less  than  those 
predicted  by  a  normal  shock.  To  empirically  account 
for  these  effects,  the  computed  normal  shock  loss  is 
reduced  by  dividing  by  the  average  inlet  relative 
Mach  number  squared. 

All  the  other  blade  row  losses— profile,  secondary, 
etc. — are  represented  by  a  correlation  with  fraction 
of  passage  height  and  aerodynamic  blade  loading. 
The  values  for  such  a  correlation  are  input  in  tabular 
form.  The  aerodynamic  blade  loading  parameter  in 
the  table  is  the  diffusion  factor  of  reference  1.  In 
equation  form  it  is 


MrVe) 
a(ri  +r2)V{ 


(2) 


The  loss  parameter  in  the  table  is 

UCOS  02 
2a 


(3) 


where  w  is  the  loss  coefficient. 

u  =  p2i  -Pi  (4) 

P\  ~Pi 

The  rotor  exit  tangential  velocity  is  calculated 
directly  from  the  Euler  equation 

(5) 

J  M 

Note  that  the  enthalpy  change  is  evaluated  by  using 
an  integral  for  the  calorically  nonperfect  gas;  that  is, 
Cp  is  a  function  of  temperature.  All  state  processes  in 
tne  program  use  thermally  perfect,  but  calorically 
nonperfect,  gas  relations;  so  integrations  and  in  some 
cases  iterations  are  used  in  several  small  function 
routines. 

Inner  loop. — The  basis  function  of  the  inner  loop 
is  to  determine  the  axial  velocity  profile  at  the 
calculation  station.  The  axial  velocity  level  is  set  by 
flow  continuity,  and  the  distribution  is  controlled  by 
the  radial  equation  of  motion.  The  differential 
equation  is  developed  in  appendix  C.  The  form  used 
in  the  program  is 
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A  velocity  gradient  procedure  is  used  to  construct 
the  axial  velocity  profile  from  the  tip  to  the  hub  with 
the  stagnation  state  values,  the  streamline 
characteristics,  and  the  tangential  component  of 
velocity  held  fixed.  Since  this  inner  loop  of  the 
program  is  used  many  times,  some  effort  was  made 
to  evaluate  its  accuracy  and  efficiency  for  typical 
streamline  spacing.  Reasonably  good  accuracy  and 
stability  were  found  to  result  from  a  rather  simple 
procedure.  Let 
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With  a  and  b  constants  for  the  /  interval  along  the 
station  path,  the  solution  for  Vm  is 


d\nP  d(rVe ) 

+  R  dl  9  rdl 


V2  —  {  -  4-  V2 
mj+  1 


') 


,2<K/-/o) 


a 

b 


(ID 


7 


t 


>#■  * 


A  two-step  procedure  is  used  in  the  program.  First 
a,  b,  and  Vm  values  on  the  streamline  j  are  used  to 
determine  a  temporary  ^mj + 1 •  The  a  and  b  values 
are  slightly  dependent  on  Vm  so  Vmj  + 1  is  used  to 
determine  new  a  and  b  values.  The  second  step  uses 
the  average  of  the  old  and  new  respective  values  of  a 
and  b  to  compute  a  final  Vmj  + 1  value.  This  1 
value  will  then  be  used  as  the  current  Vmj  value  for 
the  next  /  interval. 

When  Vm  values  are  set  on  all  streamlines,  flow 
continuity  is  checked  by  using  dr  integration  of  a 
piecewise  cubic  curve  fit  of  pVmr  values  at  the 
streamlines.  If  the  integrated  weight  flow  is  not 
within  0.01  percent  of  its  specified  value,  the  tip 
reference  Vm  is  adjusted  and  the  Vm  profile  is 
reconstructed.  The  method  of  adjusting  the  reference 
value  of  Vm  is  shown  graphically  in  figure  5.  There 
are  two  solutions  to  the  continuity  equation  in 
compressible  flow — the  subsonic  and  supersonic 
solutions.  When  a  parabolic  fit  of  trial  solutions  is 
used  to  get  a  new  trial  value  of  Vm ,  the  lower  or 
subsonic  solution  is  always  sought.  The  Vm 
adjustment  between  iterations  usualiy  is  small;  so 
convergence  normally  is  achieved  in  three  or  four 
passes. 

Once  convergence  is  achieved,  the  profile  is  back 
integrated  to  find  the  fraction  of  weight  flow  points 
represented  by  the  streamlines.  These  points  are 
saved  until  the  outer  loop  pass  through  all  the 
stations  is  completed  for  the  purpose  of  relocating 
streamlines. 

Blade  Design 

A  blade  is  defined  from  stacked  blade  elements. 
The  procedure  for  laying  out  blade  elements  and 
stacking  them  for  blade  definition  is  given  in  detail  in 


figure  5.  -  Meridional  velocity  adjustment  tor  (low  continuity 
iteration. 
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reference  2.  Only  a  summary  description  is  given 
herein.  A  blade  element  is  laid  out  on  a  cone  with  a 
center  axis  coincident  with  the  turbomachine  axis  of 
rotation.  The  angle  and  location  of  the  cone  are  fixed 
by  the  intersection  of  the  streamline  with  the  leading- 
and  trailing-edge  station  lines  of  the  blade  (fig.  6). 

The  leading-  and  trailing-edge  blade  angles  are 
related  to  aerodynamic  flow  angles  primarily  through 
two  key  correlation  parameters— incidence  angle  and 
deviation  angle.  The  user  has  some  options  for  the 
specification  of  these  correlation  parameters,  as 
already  discussed  in  the  section  on  data  input. 
Application  of  incidence  and  deviation  angles  to  the 
flow  angles  at  the  blade  edges  gives  blade  angles  in 
the  local  stream  wise  direction.  Corrections  to 
“cascade”  deviation  angle  for  a  change  in  radius  and 
axial  velocity  are  made  internally  to  the  program. 
These  corrections  are  presented  in  reference  4  to 
relate  deviation  angle  to  a  cascade  section  with 
equivalent  circulation  rather  than  with  the  same 
camber  angle. 

Because  the  cone  angle  of  the  associated  blade 
element  is  usually  a  little  different  from  these  local 
streamwise  blade  angles,  corrections  are  made  with 
current  streamwise  and  radial  direction  derivatives. 
The  blade  element  leading-  and  trailing-edge  angles 
are  calculated  from  aerodynamic  flow  angles  in 
subroutine  BLADE. 

Blade  element  layout. — There  are  several  options 
for  controlling  the  blade  element  layout  (see  the 
1DEF  (1ROW)  parameter  description  in  appendix  B). 
With  all  but  one  of  these  options  a  blade  element  is 
described  by  a  prescribed  thickness  applied  to  a 
prescribed  centerline  (fig.  4).  The  centerline  is  treated 
as  two  segments  that  are  joined  at  the  reference 
transition  point.  The  rate  of  change  of  the  local  blade 
angle  with  path  distance,  k=/(s)  (fig.  4),  is  controlled 
by  a  fourth-degree  polynomial  for  each  segment.  The 
coefficients  for  the  polynomials  are  input,  but  they 
are  scaled  in  the  program  to  match  blade  element 
inlet  and  outlet  angles.  The  fourth-degree  polynomial 
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figure  b.  I  onical  coordinate  svsleo 
tor  blade  element  layout. 


representation  of  segment  blade  angle  represents 
greater  specification  freedom  than  does  the  linear 
specification  of  reference  2,  where  the  ratio  of  inlet- 
to-outlet  segment  curvature  at  the  transition  point  is 
input  rather  than  any  polynomial  coefficients.  A 
summary  derivation  of  the  equations  for  the 
centerline  coordinates  is  given  in  appendix  D. 

Blade  element  thickness  is  defined  along  a  path 
that  is  locally  normal  to  the  centerline.  The  pressure 
and  suction  surfaces  are  equidistant  from  the 
centerline.  Thickness  is  specified  in  both  the  forward 
and  rearward  directions  from  the  maximum 
thickness  point  by  polynomials  of  the  form 
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The  input  coefficients  are  scaled  to  meet  the  leading- 
and  trailing-edge  ellipses  at  the  appropriate  tangency 
points.  The  control  routine  for  the  blade  element 
layout  in  the  program  is  CONIC. 

Blade  element  stacking.—  The  rotating  parts  of 
turbomachinery  normally  operate  at  high  stress  levels 
because  of  high  centrifugal  force.  The  high 
centrifugal  acceleration  also  causes  stress  from 
bending  moments  to  be  very  sensitive  to  blade 
element  location.  Thus  it  behooves  the  designer, 
first,  to  be  reasonably  accurate  in  the  stacking 
computation  and,  second,  to  try  to  minimize  stresses 
that  can  be  easily  reduced — namely,  those  from  the 
steady-state  bending.  The  blade  bending  moments 
from  aerodynamic  forces  can  be  counterbalanced  by 
centrifugal  force  moments  with  slight  blade  lean  in 
both  the  (r-z)  and  (r-6)  planes. 

The  reference  line  for  stacking  purposes  is  a  radial 
line  through  the  hub  stacking  reference  point  (fig.  7). 
The  sections  used  for  stacking  alignment  are  planes 
normal  to  this  reference  line  in  space.  Such  planes  are 
used  because  their  centers  of  area  are  essentially  the 
centers  of  centrifugal  force  also.  The  stacking  line  is 
a  line  that  can  be  leaned  from  the  reference  line  at  the 
hub  reference  point.  For  alignment  purposes  the 
planes  pass  through  the  stacking  line  intersection  of 
blade  elements  (fig.  7).  Blade  sections  are  defined  by 
interpolation  across  blade  elements.  When  the 
section  center  of  area  does  not  match  the  stacking 
line,  the  corresponding  blade  element  is  translated 
and  rotated  on  its  cone  for  the  stacking  adjustment. 
Normally  the  adjustments  decrease  by  about  an  order 
of  magnitude  for  successive  passes  through  the 
stacking  procedure.  For  each  pass  the  stacking  axis 
lean  angles  in  both  the  (r-z)  and  (r-9)  planes  are 


figure  f.  ■  loiiiticn  of  Blade  sections  tor  Made  element 
stacking  ailiustneuts. 


recomputed  and  adjusted  if  the  stacking  axis  lean 
option  is  activated  through  the  input  data. 

Terminal  Calculations  and  Output 

The  program  output  of  an  example  two-stage 
compressor  is  shown  in  table  III.  In  general  the 
output  is  printed  shortly  after  its  computation  so  that 
large  arrays  of  data  are  not  stored.  Data  are  printed 
from  each  of  the  major  phases  of 
computation — input,  iteration,  and  terminal.  The 
first  information  (table  III  (a))  is  the  input  data, 
which  are  printed  directly  from  input  routines  in  very 
nearly  the  order  in  which  the  input  was  read. 

The  second  major  part  of  output  (table  III  (b)>, 
from  the  iterative  phase  of  computation,  is  printed  to 
help  the  user  monitor  the  solution.  Although  these 
data  have  little  value  once  the  solution  is  converged, 
they  are  quite  helpful  in  disclosing  bad  input  and  in 
finding  sources  of  problems  when  solutions  are  not 
achieved. 

For  computational  stability  a  station  aspect  ratio, 
defined  as  (r, -r/,)/(z/+  i  -Z/),  is  limited  to  7  for 
streamline  fits.  When  the  limit  is  exceeded,  particular 
stations  (according  to  the  priorities  set  forth  in  the 
section  User  Information)  are  eliminated  from  the 
curve  fits  used  to  locate  streamlines.  The  first  data 
shown  from  the  iteration  phase  are  a  table  of  such 
calculation  station  information  (table  II  (b)).  On  the 


left  is  a  list  of  calculation  station  locations  used  to 
compute  streamlines,  along  with  the  associated 
aspect  ratios.  On  the  right  is  the  input  list  of  station 
locations  and  aspect  ratios.  When  blade  rows  are 
stacked,  the  blade  edge  stations  are  relocated,  and 
thus  the  station  aspect  ratios  change.  After  the  first 
stacking  on  iteration  4,  the  station  aspect  ratios  are 
rechecked  and  changes  in  the  station  list  are  made  if 
necessary. 

Arrays  of  axial  velocities  throughout  the  flow  field 
for  each  iteration  are  the  bulk  of  the  output  printed 
from  the  iterative  phase  of  computation.  These  data 
are  useful  for  observing  solution  stability  since  the 
solution  convergence  criterion  is  based  on  changes  of 
axial  velocity  between  successive  iterations.  Some 
compressor  overall  parameters  are  shown  above  the 
velocity  arrays.  Parameters  included  are  the  overall 
values  of  input  pressure  (PR),  current  computed 
pressure  ratio  (CPR),  enthalpy  increase  (DHC),  and 
ideal  enthalpy  increase  (DHI). 

When  the  aerodynamic  solution  is  converged,  the 
overall  parameters  for  individual  bla'de  rows  and  the 
overall  cumulative  values  in  the  compressor  are 
computed  and  printed.  Overall  temperature  and 
pressure  values  are  calculated  by  mass  averaging  their 
equivalent  enthalpy  values.  The  cumulative  forward 
axial  thrust  is  the  axial  force  exerted  on  the  rotating 
shaft  by  aerodynamic  forces  from  the  hub  inlet 
station  of  the  first  blade  row  to  the  local  point.  The 
thrust  force  shown  for  individual  blade  rows  is  the 
axial  force  on  the  shaft  from  the  trailing  edge  of  the 
upstream  blade  to  the  trailing  edge  of  this  blade  row. 
Since  the  blade  forces  on  stationary  blade  rows  act  on 
the  casing,  the  thrust  value  on  the  rotating  shaft  is 
simply  the  static  pressure  force  on  the  tapered  shaft 
in  the  forward  axial  direction.  Effects  of  cavities 
below  the  hub  flow  path  are  not  included  since 
undetermined  information  about  seal  locations  and 
pressure  differences  would  be  needed.  The  gas 
bending  moments  are  values  for  a  single  blade.  The 
bending  moments  are  referenced  to  the  stacking  axis 
intersection  with  the  flow  path  wall  from  which  the 
blades  are  attached. 

Sets  of  calculation  station  data  for  streamlines 
across  the  channel  follow  the  overall  data.  For  all 
stations,  velocity  components,  streamline  slope  and 
curvature,  and  both  stagnation  and  static  values  of 
temperature  and  pressure  are  given.  For  stations  at 
blade  row  edges,  additional  information  is  computed 
and  printed.  These  parameters  are  (1)  a  complete 
description  of  velocity  triangles,  (2)  definition  of 
blade  elements,  (3)  relations  between  aerodynamic 
and  blade  angles,  (4)  aerodynamic  performance 
parameters,  (5)  streamline  choke  area  margin,  (6) 
local  blade  force  intensity  in  pounds  per  radial  inch 
on  a  blade,  and  (7)  blade  edge  direction  derivatives 
r  dd/dr. 
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Figure  8.  -  Coordinate  system  for  made  section 
output  data. 


If  the  input  options  call  for  fabrication 
coordinates,  they  are  printed  after  all  the 
aerodynamic  output.  The  coordinates  are  printed  in 
tabular  form  with  four  sections  on  a  page,  as  shown 
in  table  III(c).  The  length  coordinate  L  is  a  distance 
along  the  chord  line,  with  the  most  forward  point 
being  zero  (fig.  8).  The  pressure-  and  suction-surface 
height  values  Hp  and  Hs,  respectively,  are  referenced 
from  the  chord  line.  Surface  height  values  are  given 
for  at  least  20  round-value  increments  of  L;  also 
surface  coordinates  are  given  for  three  specific  values 
of  L— the  blade  trailing  edge  and  the  leading-  and 
trailing-edge  ellipse  tangency  points  with  the 
surfaces. 

A  blade  section’s  properties  are  shown  above  its 
table  of  coordinates  (table  111(c)).  The  blade  section 
radial  location,  the  L  and  H  stacking  point  values, 
and  the  section  setting  angle  are  given  to  locate  and 
orient  the  blade  section.  The  blade  section  center-of- 
area  coordinates,  section  area,  minimum  and 
maximum  moments  of  inertia  through  the  center 
area,  orientation  angle  of  the  maximum  moment  of 
inertia  with  respect  to  the  axial  direction,  section 
torsion  constant,  and  twist  stiffness  are  all  useful 
information  for  design  and  stress  analysis. 

After  all  the  fabrication  coordinates  for  a  given 
blade  row  are  printed,  the  blade  section  coordinates 
are  presented  in  another  orientation  that  may  be 
more  useful  for  further  flow  analysis.  With  a 
stacking  axis  reference,  coordinates  for  the  same 
blade  sections  are  given  in  the  axial  and  tangential 
directions. 


User  Information 

Since  earlier  sections  of  the  report  discuss  the 
input,  output,  and  main  centers  of  program  control, 
this  discussion  is  directed  at  the  user  who  is  trying  to 
get  the  program  on  his  computer  and  to  make  it  run 
efficiently.  Some  facts  about  the  program  as  well  as 


some  advice  about  the  input  are  given. 

The  code,  which  is  written  in  FORTRAN,  takes 
about  80,000  decimal  words  of  computer  storage. 
The  call  relation  among  the  subroutines  is  shown 
graphically  in  figure  9.  Note  that  the  tickmarks  on 
the  routine  boxes  in  the  figure  mean  that  there  are 
other  call  lines  to  the  routine.  These  lines  are  shown 
on  the  other  part  of  the  figure  where  the  routine 
name  is  repeated.  The  program  running  time  on 
either  a  Univac  1110  or  an  IBM  360-67  is  about  2 
minutes  for  a  single-stage  compressor  and  about  5 
minutes  for  a  five-stage  compressor.  Several  of  the 
key  indices  in  COMMON/SCALAR/  are  described 
in  the  following  tabulation. 


Index  Description 

1  calculation  station  index  after 

preliminary  calculations  are 
completed.  The  program  is 
dimensioned  for  50  calculation  stations 
and  20  blade  rows,  of  which  only  10 
can  be  rotors.  Each  blade  row  accounts 
for  two  calculation  stations — one  at 
the  leading  edge  of  the  blade  and  the 
other  at  the  trailing  edge.  Rotors, 
stators,  and  annular  calculation 
stations  can  be  put  together  in  any 
combination  with  the  following 
constraints:  The  number  of  stations 
cannot  exceed  50.  There  must  be  at 
least  four  annular  stations  ahead  of  the 
first  blade  row  and  at  least  three 
annular  stations  behind  the  last  blade 
row. 

IROTOR  rotor  index 


I  ROW  blade  row  index 


J 

K 


streamline  index.  Streamlines  are 
numbered  from  one  at  the  tip. 

loss  set  index  for  subroutine  INPUT 


As  indicated  in  the  table  at  least  four  annular  stations 
are  expected  upstream  of  the  first  blade  row  and  at 
least  three  downstream  of  the  last  blade  row. 
Additional  annular  stations  can  be  located  between 
blade  rows  but  not  within  blade  rows;  that  is,  not 
between  the  inlet  and  outlet  stations  of  a  given  blade 
row. 

Streamline  intersections  of  station  lines  are 
determined  by  integrating  velocity  profiles  at  station 
lines  to  the  specified  mass  flow  fractions.  Streamline 
slope  and  curvature  are  determined  from  streamwise 


curve  fits  of  these  intersections.  The  consequence  of 
this  procedure  is  that  the  number  of  iterations  and 
the  program  convergence  characteristics  are 
dependent  on  the  calculation  station  location 
although  the  final  solution,  in  general,  is  not  very 
dependent  on  the  location  of  the  calculation  stations. 

The  user  can  reduce  the  number  of  iterations  and 
hence  the  program  running  time  with  good  placement 
of  calculation  stations.  The  first  calculation  station 
should  be  placed  upstream  of  the  first  blade  row  a 
distance  at  least  equal  to  two  or  three  annulus 
heights.  The  best  far-upstream  inlet  condition  is 
straight  axial  flow  with  no  wall  curvature.  Less 
iterations  are  usually  needed  for  more  widely  spaced 
calculation  stations;  however,  enough  iterations 
should  be  used  to  properly  locate  the  streamlines. 
Calculation  station  spacing  can  vary  somewhat  along 
the  annulus  but,  as  a  general  guideline,  successive 
station  increments  should  not  be  changed  more  than 
35  percent. 

When  calculation  stations  are  input  close  together, 
only  some  of  them  will  be  used  for  locating  the 
streamlines  if  the  station  aspect  ratio  is  above  7.0. 
This  is  done  for  program  stability  and  convergence 
toward  a  solution.  If  the  user  does  not  specify  which 
stations  to  eliminate  from  the  streamline  location 
procedure,  the  program  has  logic  to  do  so  when  the 
station  aspect  ratio  exceeds  7.0.  The  priority  of 
stations  kept  for  streamline  location  is  as  follows:  (1) 
blade  row  exit  stations  are  always  used,  (2)  blade  row 
inlet  stations  are  kept  if  the  blade  row  aspect  ratio  is 
less  than  7.0,  and  (3)  an  annular  station  is  kept  if 
neither  adjacent  station  is  closer  than  the  aspect  ratio 
tolerance. 

The  user  can  also  specify  that  particular  annular 
stations  not  be  used  for  streamline  definition  through 
the  alphanumeric  station  designation.  The  program 
looks  for  ROTO  for  rotor,  STAT  for  stator,  or 
ANNU  for  regular  annular.  Any  other  combinations 
of  letters,  numbers,  or  symbols  designates  the  station 
as  the  extra-annular  type.  All  the  computations  that 
are  done  for  regular  annular  stations  also  are  done 
for  the  extra-annular  stations.  The  only  difference  is 
that  the  new  streamline  locations  at  that  station  are 
not  used  for  the  curve  fit  for  streamline  parameters. 
When  the  new  curve  fit  streamlines  are  established, 
their  intersections  with  the  station  line  are  found  and 
the  streamline  parameters  at  that  point  are  used  in 
the  equation-of-motion  calculations. 

The  arrays  of  points  that  describe  the  hub  and  tip 
casing  contours  should  extend  at  least  from  the 
furthest  upstream  calculation  station  to  the  furthest 
downstream  one.  There  should  be  enough  data 
points  to  adequately  define  the  desired  casing 
contours  with  a  spline  curve  fit. 

The  input  boundary  layer  blockage  factors  have  an 
option.  A  displacement  thickness  from  the  wall  can 


it 


be  specified  instead  of  blockage  as  a  fraction  of 
annulus  height.  This  is  done  by  using  a  negative 
number  the  magnitude  of  which  is  the  value  of 
displacement  thickness. 

A  total  pressure  profile  can  be  input  in  place  of 
losses.  Although  the  way  to  activate  this  option  has 
been  discussed  earlier,  its  full  effects  need  to  be 
understood.  This  option  is  activated  for  a  particular 
blade  row  by  using  zero  or  a  negative  number  in 
ILOSS  (1ROW).  When  the  option  is  activated,  an 
additional  data  card  is  required  for  that  blade  row 
(fig.  12(a)).  The  first  parameter  PTT(IROW),  or  P, 
in  the  equation,  is  the  blade  row  tip  (larger  radius) 
total  pressure  in  psia.  The  five  other  parameters  are 
polynomial  constants  P\  to  P$\  therefore  a  total 
pressure  at  some  other  radial  location  is 

P  =  P,(1.0  +  P,P  +  P2R2  +P3P3  +  P4/?4  +  P5P5) 


where 


r,-rh 

or  the  fraction  of  passage  height  at  the  blade  row 
exit.  Because  these  coefficients  are  stored  into  the 
locations  of  loss  sets  4  and  5,  those  loss  sets  are 
destroyed  for  the  run  even  if  read  in. 

When  the  pressure  level  is  specified  instead  of 
losses  for  the  last  blade  row  of  the  compressor,  there 
is  an  overspecification  of  data  because  the  inlet 
pressure  and  compressor  pressure  ratio  are  input  too. 
In  computation  the  pressure  ratio  predominates;  so 
the  pressure  levels  will  be  adjusted  as  necessary.  Also 
note  that  when  the  pressure  level  is  input,  the  total 
temperature  profile  must  also  be  input  (table  II). 
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lal  Subroutines  used  in  input  and  iteration  phases. 
Figure  9.  -  Line  representation  of  subroutine  calls. 


Figure  9.  -  Concluded. 


At  a  rotor  exit  the  total  temperature  level  can  be 
input  in  place  of  the  cumulative  energy  addition 
fraction.  If  the  input  CRENGY  (IROTOR)  is  greater 
than  2.0,  the  value  is  interpeted  to  be  the  rotor  exit 
tip  temperature  in  degrees  Rankine.  In  the 
preexecution  phase  of  computation  the  temperature 
is  converted  and  used  as  an  appropriate  energy 
addition  value.  The  polynomial  coefficients  for  the 
radial  distribution  of  total  temperature  are  input  in 
the  former  pressure  polynomial  coefficient  locations, 
PARA(IROW)...PRE(lROW).  During  regular 
iteration  the  program  will  use  the  polynomial  form 


for  rotor  exit  total  temperature  distribution  when 
|PRA(IROW)]>  100.0.  The  polynomial  coefficient 
represented  by  PRA(IROW)  is  found  by  adding  or 
subtracting  the  number  of  100’s  needed  to  give  a 
remainder  in  the  range  -  100.0  to  100.0. 

When  the  total  temperature  level  is  input,  the  total 
pressure  level  can  be  set  in  two  ways.  It  can  either  be 
determined  from  losses  or  input  directly  by  a 
polynomial,  as  discussed  earlier  in  this  section. 

The  description  of  parameter  variations  with 
polynomials  assures  smoothness,  but  the 
specification  of  polynomial  coefficients  is  not  always 
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easy.  In  most  cases  the  range  of  applicability  for  the 
polynomial  independent  variable  is  0  to  1.0.  This 
considerably  eases  the  burden  on  the  user  since 
computation  is  normally  not  needed  to  choose  and 
set  the  polynomial  coefficients.  When  the  higher 
degree  terms  are  used  to  define  distributions,  the  end 
conditions  are  relatively  easy  to  meet.  However, 
some  simple  computations  are  needed  to  check  the 
distribution. 

Another  caution  is  that  combinations  of 
reasonable-looking  numbers  often  give  blade 
elements  that  one  can  judge  to  be  poor  by  visual 
observation.  The  capability  to  make  machine  graphic 
plots  of  blade  elements  and  the  channel  formed  by 
adjacent  blades  is  very  useful.  Such  plots  are  made  in 
subroutine  EPLOT,  which  is  activated  by  the  input 
parameter  OPM.  Since  graphics  packages  differ  with 
computer  systems,  the  program  presented  will  not 
necessarily  work  directly  on  a  user's  computer. 
However,  it  is  suggested  that  the  user  make  the 
conversions  necessary  to  plot  the  blade  element 
surface  arrays  generated  in  EPLOT. 

The  determination  of  acceptable  polynomial 
coefficients  for  the  centerline  and  thickness  of  an 
entire  blade  row  can  be  difficult  when  high-degree 
terms  are  used.  This  task  was  eased  considerably  at 
NASA  Lewis  with  an  interactive  graphics  capability. 
A  series  of  computer  programs  were  developed  to 
design  particular  blade  elements  from  actual 
centerline  angle  and  thickness  distributions.  These 
data  were  then  curve  fit  by  least-squares  methods  to 
produce  the  input  required  by  the  program  described 
in  this  report.  Visual  observaton  of  blade  elements 
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generated  by  this  input  for  several  fractions  of 
annulus  height  is  very  helpful  in  avoiding  obviously 
unacceptable  configurations. 

The  computer  peripheral  equipment  also  can  be 
used  by  some  other  subroutines  when  options  are 
activated  with  the  parameter  OPO.  When  the  punch 
option  is  activated,  the  tables  of  fabrication 
coordinates  shown  on  the  listing  are  punched  on 
cards  in  subroutine  COORD.  When  the  plot  option 
of  OPO  is  activated,  subroutine  BLUEPT  plots 
tables  of  fabrication  coordinates  on  a  blueprint 
format.  If  a  plot  option  is  activated  by  either  OPM  or 
OPO,  subroutine  MERID  is  also  called.  It  produces  a 
meridional  plane  plot  of  the  annulus  flow  path  with 
the  calculation  stations  and  streamlines  included. 

This  code  is  interfaced  with  three  other  NASA 
codes  through  punched  card  output.  Input  for  the 
TSONIC  code  (ref.  5),  which  is  a  blade-to-blade 
channel  flow  analysis  code,  is  obtained  with  the  T 
option  of  OPM.  Input  for  the  MERIDL  code  (ref.  6), 
which  is  a  more  detailed  hub-to-shroud  flow  analysis 
code  within  a  blade  row,  is  obtained  with  the  M 
option  of  OPO.  Input  for  an  off-design  performance 
prediction  code  that  is  being  developed  at  NASA 
Lewis  is  obtained  with  the  O  option  of  OPM. 

The  computer  program  can  be  obtained  from 
COSMIC,  112  Barrow  Hall,  University  of  Georgia, 
30601.  The  COSMIC  program  number  is 
LEW-13505. 

Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  December  29,  1980 


Appendix  A 
Symbols 


A  annulus  area;  also  streamtube  channel  area 

Aj  polynomial  constants  for  as  a  function  of  S 
a  sonic  velocity,  ft/sec;  also  a  coefficient  in 

velocity  gradient  equation;  also  a 
polynomial  coefficient 

b  coefficient  in  velocity  gradient  equation; 

also  a  polyf  -mial  coefficient 
C  constant 

C,  polynomial  constants  for  conic  radius  as  a 
function  of  S 

c  blade  chord,  in.;  also  a  polynomial 

coefficient 

cp(t)  specific  heat  function  for  constant 

pressure,  ft/sec2  °R 
D  blade  element  diffusion  factor 

Dj'j- |i00  simplified  nomenclature,  D,  =  - (Ct)/(i)R, 
d  polynomial  coefficient 

/  friction  force,  ft/sec2 

H  stagnation  enthalpy,  ft/sec2 

Hp  pressure-surface  height,  in. 

Hs  suction-surface  height,  in. 

h  static  enthalpy,  ft/sec2 

/  integer  index;  also  incidence  angle,  deg 

j  integer  index 

k  curvature  in  curvilinear  coordinate  system, 

ft '  1 ;  also  an  integer  index 
L  distance  along  chord  line,  in. 

/  distance  along  calculation  station  line,  in. 

M  Mach  number 

m  streamline  direction  in  meridional  plane, 

in.;  also  an  integer  index 

n  streamline  normal  direction  in  meridional 

plane,  in. 

P  stagnation  pressure,  lb/ft2 

p  static  pressure,  lb/ft2 

R  conic  coordinate  radius,  in. 

Rn- 1  oo  series  coefficients  for  polynomial, 
R,/R  =  l  +R{S  +  R2S2  +  R3S3  +  ... 

Rm  radius  of  curvature  in  meridional  plane,  ft 
(R  gas  constant,  ft  lb/slug  °R 

r  radius  from  axis  of  rotation,  in. 

S  blade  element  path  distance,  in. 

s  entropy,  ft/sec2  °R 

T  stagnation  temperature,  "R 

t  static  temperature,  °R;  also  blade  element 

thickness,  in. 


U  local  blade  velocity,  ft/sec 

u  generalized  variable  in  a  differential 

equation 

V  velocity,  ft/sec 

v  generalized  variable  in  a  differential 

equation 

weight  flow,  lb/sec 
z  axial  distance,  in. 

a  angle  of  streamline  with  reference  to  axial 

direction,  deg 

13  flow  angle  relative  to  meridional  direction, 

deg 

7  blade  chord  angle,  deg 

6  deviation  angle,  deg 

f  angular  coordinate  on  blade  element  layout 

cone,  rad 

6  circumferential  direction,  rad 

k  blade  angle  relative  to  local  conic  ray,  deg 

\  local  angle  of  calculation  station  line  with 

reference  to  radial  direction,  deg 
p  static  density,  slug/ft3 

a  blade  element  solidity,  chord/tangential 

spacing 
r  time,  sec 

w  loss  coefficient 

Subscripts: 

ca  center  of  area 

/  calculation  station  index 

i  ideal  value,  as  by  an  isentropic  process 

j  streamline  index 

le  leading  edge 

m  streamline  direction  in  meridional  plane; 

also  maximum  thickness 

n  streamline  normal  direction  in  meridional 

plane 

o  initial  value 

sp  stacking  point 

t  transition  point 

te  trailing  edge 

d  circumferential  direction 

1  blade  row  inlet 

2  blade  row  outlet 

Superscript: 

( )'  relative  to  rotor 

(  )*  flow  at  sonic  condition  (Af'  =  1.0) 


Appendix  B 

Input  Parameters  for  Compressor  Design  Program 

The  input  variables  for  the  compressor  design  program  and  the  associated  options  are  described  in  this 
appendix.  The  format  for  the  input  data  is  given  in  figures  10  to  12.  The  calculation  station  and  blade  row  data 
sets  are  input  in  the  order  in  which  they  occur  in  the  compressor  flow.  If  any  of  the  sets  of  option  cards  for 
blade  rows  are  needed,  they  are  considered  part  of  the  blade  row  set  and  they  follow  the  particular  basic  blade 
row  data  set  in  the  order  shown  in  figure  12.  The  only  exception  is  any  XCUT  cards  that  are  read  in  the  output 
routines.  These  cards  are  at  the  end  of  the  input  data,  but  of  course  the  sets  of  XCUT  values  must  be  placed  in 
the  same  order  as  the  stations  specifying  them. 

In  the  following  list  of  parameters  the  independent  variable  S  appears  frequently.  Since  it  is  an  important 
blade  element  definition  variable,  this  preliminary  explanation  of  its  definition  and  usage  is  given.  The  variable 
S  in  equations  for  the  blade  element  centerline  is  the  distance  in  either  direction  from  the  transition  point  as  a 
reference.  The  variable  S  in  equations  for  the  thickness  distribution  is  the  distance  in  either  direction  from  the 
maximum  thickness  point  as  a  reference.  All  four  of  these  usages  of  S  are  shown  in  figure  4.  In  all  cases,  S 
values  are  positive  away  from  their  reference  point.  The  S  values  for  thickness  definition  are  normalized  by 
blade  element  chord.  The  S  values  for  centerline  definition  are  also  normalized  by  blade  element  chord  when 
IDEF(IROW)  is  less  than  zero;  however,  when  IDEF(IROW)  is  greater  than  zero,  S  is  normalized  to  1.0;  that 
is,  the  maximum  segment  S  is  1.0. 


Figure  10.  -  Input  date  format  ol  qeneral  information. 
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Figure  11.  -  Input  data  format  of  calculation  stations  and  basic  blade  row  information. 


Parameter 


Description 


Formal 


AA 


AB 


This  parameter  is  used  twice  to  indicate  options  in  alphanumeric  form.  As  the  A4 
first  term  of  a  data  set  it  indicates  the  type  of  calculation  station  or  blade  row 
(ANNULAR,  ROTOR,  or  STATOR).  Any  station  description  other  than 
ANNU,  ROTO,  or  ST  AT  will  be  treated  as  an  extra-annular  station,  that  is,  the 
streamlines  will  not  be  forced  to  pass  through  the  streamtube-fraction-of 
weight-flow  point  as  determined  by  continuity  at  the  station.  The  second  use  of 
AA  later  in  the  data  set  is  the  incidence  angle  option  for  blade  design  purposes. 
Interpretable  options  are  2-D,  3-D,  SUCTION,  and  TABLE.  A 
noninterpretable  incidence  option  word  is  set  to  the  2-D  option.  The  2-D  and 
3-D  options  mean  incidence  angles  are  determined  by  procedures  in  reference  1 
for  the  respective  option.  The  suction  option  gives  zero  incidence  to  the  suction 
surface  of  the  blade  at  the  leading  edge.  The  TABLE  option  means  the  blade 
incidence  angles  for  the  blade  element  will  be  input  in  tabular  form, 
INC(IROW,J),  at  the  end  of  the  data  set. 

This  parameter  completes  the  incidence  TABLE  option  discussed  above.  To  A4 
reference  incidence  lo  the  suction  surface  at  the  leading  edge,  the  eight  spaces  of 
the  card  for  AA  and  AB  must  read 

TABLE  SS 
AA  AB 


ACF(l.IROW), 

ACF(2,lROW), 

ACF(3,IROW), 

ACF(4,IROW) 


(If  AB  is  anything  other  than  E  SS,  the  incidence  angles  will  be  referenced  to  the 
leading-edge  centerline.) 

polynomial  coefficients  for  linear  coefficient  of  blade  element  centerline  angle  F10.4 
equation  for  front  segment,  k=k,+oS  +  bS2  +cS2  +dS*  with 
o  =  ACF1  +  ACF2*R  +  ACF3*R2  +  ACF4•R-^  where  R  =  (r, -r)/(r, -re¬ 
fraction  of  passage  height  at  blade  leading  edge 


ACR(l,IROW),  same  as  above  for  rear  segment  with  same  R 
ACR(2,IROW), 

ACR(3,IROW), 

ACR(4,IROW) 


F10.4 
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Figure  12.  -  Input  data  format  of  additional  bfade  row  information  if  needed  by  the  options. 


Parameter  Description  Format 

ALIM(IROW)  For  a  data  set  designated  ROTOR,  ALIM(IROW)  is  the  minimum  allowable 
relative  flow  angle  (deg)  leaving  the  rotor  hub.  For  a  data  set  designated 
STATOR,  ALIM(IROW)  is  the  maximum  Mach  number  entering  the  stator  at 
the  hub.  The  program  will  reduce  the  stage  energy  addition  to  satisfy  these 
conditions  if  a  limit  criterion  has  been  reached  during  computation.  If  no 
aerodynamic  limits  have  been  reached  in  some  other  stages  of  a  multistage 
compressor,  the  program  will  try  to  pick  up  the  energy  loss  of  the  limiting  stage 
in  the  stages  free  of  aerodynamic  limits.  If  all  stages  have  reached  some 
aerodynamic  limit,  the  overall  compressor  pressure  ratio  is  degraded  to  get  all 
stages  within  the  specified  aerodynamic  limits.  The  most  efficient  way  to  run  the 
program  is  to  specify  the  stage  energy  addition  levels  so  than  aerodynamic  limits 
%  are  not  reached  or  at  least  not  reached  in  a  drastic  fashion. 

/* 
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t 


Parameter 


Description 


Format 


ATF(l.IROW), 

ATF(2,IROW), 

ATF(3,IROW), 

ATF(4,IROW) 


ATR(l.IROW), 

ATR(2,lROW), 

ATR(3,IROW), 

ATR(4,IROW) 

BB 


BCF(l.IROW), 

BCF(2,IROW), 

BCF(3,IROW), 

BCF(4,lROW) 

BCR(l,IROW), 

BCR(2,IROW), 

BCR(3,IROW), 

BCR(4,IROW) 

BH(I) 


BMATL(IROTOR) 


BLADES(IROW) 

BLEED(I) 

BT(I) 


polynomial  coefficients  for  first  coefficient  a  of  blade  element  thickness  equation  FI0.4 
forward  of  maximum  thickness  point 

i  =  Yc  ~a(^o -S- VS0  +  JL)  -  bS2  -  cS2  -  dS 4 

with  a  =  ATF1  +  ATF2*/?  +  ATF3»/?2  +  ATF4»fl3,  where  R  is  fraction  of 
passage  height  at  blade  leading  edge  and  Sa  is  distance  from  maximum 
thickness  point  to  centerline  intersection  of  edge  ellipse  (fig.  4) 

same  as  above  for  rearward  thickness  with  same  R  F10.4 


deviation  angle  option  for  blade  design  purposes.  Interpretable  options  are  2-D,  A4 
3-D,  TABLE,  CARTER,  and  MODIFY.  Noninterpretable  input  is  set  to  the 
2-D  option.  For  the  2-D  and  3-D  options,  deviation  angles  are  determined  by 
procedures  of  reference  1  for  the  corresponding  option.  The  CARTER  and 
MODIFY  options  are  now  the  same  in  the  program.  They  indicate  the  use  of  a 
Carter’s  rule  with  a  modification  when  the  front  and  rear  segments  of  a  blade 
element  have  different  camber  rates.  The  TABLE  option  means  that  the  blade 
deviation  angles  for  the  blade  elements  will  be  input  in  tabular  form, 
DEV(IROW,J),  at  the  end  of  the  data  set. 

polynomial  coefficients  for  quadratic  coefficient  of  blade  element  centerline  F10.4 
angle  equation  for  front  segment,  k  =  k  ,  +  aS  +  bS2  +  cS2  +  dS* 
with  fc  =  BCFl  +BCF2*/?  +  BCF3»/?2  +  BCF4*/?3,  where 
R  =  (r,-r)/(r,  ~rh)— fraction  of  passage  height  at  blade  leading  edge 

same  as  above  for  rear  segment  with  same  R  F10.4 


hub  blockage  factor  for  each  calculation  station;  fraction  of  the  station  annular  F10.4 
area  to  be  allowed  for  hub  annular  surface  boundary  layer  blockage.  The  hub 
streamline  will  be  displaced  away  from  the  physical  wall  a  distance  that  gives  the 
specified  annular  fraction.  Negative  input  values  are  used  as  the  magnitude  of 
boundary  layer  displacement  in  inches. 

rotor  material  density  (lb/in3).  If  a  positive  nonzero  number  is  input,  the  blade  F10.4 
will  be  stacked  so  as  to  balance  out  gas  bending  moments  with  the  centrifugal 
force  moment  for  the  material  density.  Because  the  hub  stacking  point  stays 
fixed,  the  tip  location  is  moved  if  necessary. 

number  of  blades  in  each  rotor  or  stator  blade  row  F10.4 

fraction  of  weight  flow  bled  off  at  particular  calculation  station  F10.4 

same  as  BH(I)  except  applicable  at  lip  F10.4 


Parameter 


Description 


Format 


BTF(l,lROW) 


BTR(l.IROW), 

BTR(2,IROW), 

BTR(3,IROW), 

BTR(4,IROW) 

CC 


CCF(l.IROW), 

CCF(2,IROW), 

CCF(3,IROW), 

CCF(4,IROW) 

CCR(l.IROW), 

CCR(2,IROW), 

CCR(3,IROW), 

CCR(4,IROW) 

CHORDA(IROW), 

CHORDB(IROW), 

CHORDC(IROW) 


CHOKE(IROW) 


CPCO(I) 
for  1  =  1,6 


polynomial  coefficients  for  quadratic  coefficient  of  blade  element  thickness 
equation  forward  of  maximum  thickness  point 

~  =  ^  +a(\'S^S-\IS0+  - bS2  - cS2  - dS4 

2c  2c  \  0  °  2\ISa ) 

with  b  =  BTF1  +  BTF2*/?  +  BTF3*/?2  +  BTF4*/?3,  where  R  is  fraction  of 
passage  height  at  blade  leading  edge. 

same  as  above  for  rearward  thickness  with  same  R 


blade  element  geometry  option  for  blade  design  purposes.  Interpretable  options 
are  CIRCULAR,  OPTIMUM,  and  TABLE.  The  CIRCULAR  option  gives 
circular  arc  blade  elements.  Noninterpretable  input  will  be  set  to  the 
CIRCULAR  option.  The  OPTIMUM  option  means  that  the  ratio  of  blade 
element  segment  turning  rates  will  be  set  by  an  empirical  function  of  inlet 
relative  Mach  number.  Below  an  M |  of  0.8  the  blade  element  will  be  a  circular 
arc.  As  M\  is  increased,  the  ratio  of  front  segment  turning  rate  to  rear  segment 
turning  rate  is  reduced.  A  limit  of  zero  camber  on  the  suction  surface  of  the 
front  segment  is  approached  at  an  A/j  of  about  1 .60.  The  TABLE  option  means 
the  ratio  of  blade  segment  turning  rates  will  be  input  in  tabular  form, 
PHI(IROW.J),  at  the  end  of  the  data  set. 

polynomial  coefficients  for  cubic  coefficient  of  blade  element  centerline  angle 
equation  for  front  segment,  k=k,  +  aS  +  bS2  +  cS3  +  dS4 
with  c  =  CCF1  +  CCF2 */?  +  CCF3 *R2  +  CCF4 *R\  where 
/?  =  (/-,  ~r)/(r,  —rf,) — fraction  of  passage  height  at  blade  leading  edge 

same  as  above  for  rear  segment  with  same  R 


constants  to  define  ratio  of  blade  element  chord  to  tip  chord  on  projected 
plane 

—  =  1  +  R*CHORDA(IROW)  +  R2  *CHORDB(IROW) 

ct  ip 

+  /?3*  CHORDC(IROW) 

where  /?  =  (r,  -r)/(r,  -/>)— fraction  of  annulus  height  at  blade  stacking  line 

desired  minimum  value  of  (A /A*)-  1.0,  where  A /A*  is  the  ratio  of  local 
streamtube  area  in  the  channel  to  the  area  required  when  M'  =  1.0  within  a 
blade  passage.  If  zero  is  input,  no  adjustment  will  be  attempted  within  the 
program.  For  input  values  greater  than  zero, incidence  angle  will  be  increased  as 
necessary  up  to  a  maximum  of  +  2.0°  on  the  leading  edge  of  the  suction  surface 
in  an  attempt  to  give  the  specified  choke  margin  at  the  covered  channel  entrance 
if  the  minimum  occurs  at  the  channel  inlet. 

constants  for  specific  heat  polynomial  function  of  temperature 

cp  =  CPCO(  1 )  +  CPCO(2)*  T  +  CPCO(3)»  T2  +  CPCO(4)*  P 

+  CPCO(5)*  P  +  CPCO(6)*  P 


F10.4 


F10.4 


A4 


F10.4 


F10.4 


FI  0.4 


E20.8 


Parameter 


Description 


Format 


CRENGY 

(IROTOR) 


CTF(l,IROW), 

CTF(2,IROW), 

CTF(3,IROW), 

CTF(4,IROW) 


CTR(l.IROW), 

CTR(2,IROW), 

CTR(3,IROW), 

CTR(4,IROW) 

DCF(l,IROW), 

DCF(2,IROW), 

DCF(3,IROW), 

DCF(4,IROW) 

DCR(l,IROW), 

DCR(2,IROW), 

DCR(3,IROW), 

DCR(4,IROW) 

DD 


DEV(IROW,J) 


DFTAB(K,J,I) 


desired  cumulative  energy  addition  fraction  through  particular  rotor  to  total 
energy  addition  of  compressor.  Thus  the  fractions  are  progressively  larger 
positive  numbers  through  successive  rotors.  The  last  rotor  must  have 
CRENGY  =  1.0  to  meet  the  input  pressure  ratio.  If  a  value  greater  than  2.0  is 
input,  the  value  is  interpreted  as  a  rotor  exit  total  temperature  level  in  degrees 
Rankine  instead  of  the  cumulative  energy  addition  fraction.  In  the  preexecution 
phase  of  computation  the  input  temperature  is  converted  and  used  as  an 
appropriate  energy  addition  value. 

polynomial  coefficients  for  cubic  coefficient  of  blade  element  thickness 
equation  forward  of  maximum  thickness  point 

~  =  +  a(\iS^-S-'JS0  +  -?L)-bS2 -cS*  -dS* 

2c  2c  \  0  °  2  \'S0  / 

with  c  =  CTFl +CTF2»/?  +  CTF3*/?2 +CTF4»/?3,  where  R  is  fraction  of 
passage  height  at  blade  leading  edge 

same  as  above  for  rearward  thickness  with  same  R 


polynomial  coefficients  for  fourth  degree  coefficient  of  blade  element  centerline 
angle  equation  for  front  segment,  k  =  k,  +aS  +  bS 2  +  cS3  +r/.S4 
with  d  =  DCF1  +  DCF2*/?  +  DCF3*/?2  +  DCF4*/?3,  where 
R  =  (r,  -r)/(r,-rh)— fraction  of  passage  height  at  blade  leading  edge 

same  as  above  for  rear  segment  with  same  R 


option  control  of  location  of  transition  point  between  segments  of  a  blade 
element.  The  interpretable  options  are  CIRCULAR,  SFIOCK,  and  TABLE. 
The  SHOCK  option  locates  the  transition  point  on  the  suction  surface  at  the 
normal  shock  impingement  point  from  the  leading  edge  of  the  adjacent  blade. 
The  TABLE  option  means  the  location  of  the  transition  point  will  be  input  in 
tabular  form,  TRANS  (1ROW.J),  at  the  end  of  the  data  set.  The  CIRCULAR 
option  and  noninterpretable  data  put  the  transition  point  at  midchord. 

deviation  angle  (deg)  that  can  be  specified  by  option.  If  the  tabular  option  is 
used,  a  value  is  expected  for  each  streamline  starting  from  the  tip. 

blade  element  diffusion  factor  (D  factor)  for  which  profile  losses  are  tabulated. 
Five  values  are  input  for  each  streamline;  that  is,  K  always  has  values  from  1  to 
5,  J  is  the  streamline  index,  and  1  is  the  loss  set  index.  The  maximum  number  of 
sets  is  5.  Because  D-factor  values  normally  fall  between  0.3  and  0.7,  values  of 
0.3,  0.4,  0.5,  0.6,  and  0.7  for  DFTAB  on  a  streamline  can  be  implied  by  leaving 
the  DFTAB  values  blank.  As  a  consequence  of  this  option  the  DFTAB  cannot 
be  exactly  0.0  when  K  =  1  if  you  do  not  want  the  implied  values  of  DFTAB. 

aerodynamic  D-factor  limit.  In  a  data  set  designated  ROTOR  this  limit  applies  at 
the  tip  streamline.  For  a  STATOR  data  set  the  limit  applies  at  the  hub.  The 
program  operates  with  this  limit  criterion  in  the  same  way  as  it  did  with 
ALIM(IROW). 


F10.4 

F10.4 

F10.4 

F10.4 

F10.4 

A4 

F10.4 

F8.4 

F10.4 


DLIM(IROW) 


Parameter 


Format 


Description 


DLOS(K,J,l) 


DTF(l,lROW), 

DTF(2,IROW), 

DTF(3,lROW), 

DTF(4,lROW) 


DTR(l.IROW), 

DTR(2,lROW), 

DTR(3,IROW), 

DTR(4,IROW) 

EB 


EE 


ELE(UROW), 

ELE(2,IROW), 

ELE(3,IROW), 

ELE(4,IROW) 


ETE(l.IROW), 

ETE(2,IROW), 

ETE(3,IROW), 

ETE(4,IROW) 


FLOFRA(I) 


profile  loss  parameter  <•>  cos  /ij /2a  corresponding  to  DFTAB(K,J,1)  reference 
arrays 

polynomial  coefficient  for  fourth  coefficient  of  blade  element  thickness  equation 
forward  of  maximum  thickness  point 

Tc  =  'fc  JL  ybS2  -cS3  -dS< 

with  d=DTFl  +  DTF2*/?  +  DTF3*/?2  +  DTF4*/?3,  where  R  is  fraction  of 
passage  height  at  blade  leading  edge 

same  as  above  for  rear  segment  with  same  R 


EB  completes  FABLE  option  of  maximum  thickness  location.  If  the  eight  spaces 
controlling  the  option  appear  as 

TABLE  LE 
EE  EB 

the  input  values  of  ZMAX(IROW.J)  will  be  used  as  the  fraction  of  chord 
distance  from  the  leading  edge.  If  EB  is  not  as  shown,  the  values  of 
ZMAX(IROW,J)  will  be  used  as  the  fraction  of  chord  distance  behind  the 
transition  point. 

option  control  of  location  of  maximum  thickness  point  of  a  blade  element.  The 
interpretable  options  are  TRAN  and  TABLE.  The  TRAN  option  and 
noninterpretable  options  will  set  the  maximum  thickness  point  at  the  transition 
point.  The  TABLE  option  means  the  maximum  thickness  point  location  will  be 
input  in  tabular  form,  ZMAX(IROW,J),  at  the  end  of  the  data  set. 

coefficients  for  leading-edge  ellipse 
ratio  of  semimajor  to  semiminor  axes 
minus  1 

e=--\=  ELE1  +  ELE2*fl  +  ELE3*/?2  +  ELE4*/?3 
a 

where  R  is  fraction  of  passage  height  at  blade  leading 

coefficients  for  trailing-edge  ellipse  ratio  of  semimajor  to  semiminor  axes  minus  1 

<?=  -  -  1  =ETE1  +ETE2*/?  +  ETE3*/?2  +  ETE4*R3 
a 

where  R  is  fraction  of  passage  height  at  blade  trailing  edge 

cumulative  weight-flow  split  between  streamlines  starting  from  tip.  NTUBES, 
which  is  NSTRM-1,  values  are  read.  Thus  the  first  value  is  greater  than  zero  and 
succeeding  values  must  increase  to  1 .0  in  order  for  the  last  value  to  account  for 
the  accumulation  of  flow  for  all  streamtubes. 

mass  flow  (lb/sec)  entering  the  first  calculation  station 


F8.4 

F10.4 

F10.4 

A4 


A4 

F10.4 


F10.4 


FLOW(I) 


F10.4 


Parameter 


Description 


Format 


IDEF(IROW) 


ILOSS(IROW) 


INC(IROW,J) 

MOLE 

NA 

N  BROWS 

NHUB 

NLOSS 

NTIP 

NXCUT 


NSTRM 

OP 


blade  definition  index.  When  the  index  is  zero,  the  blade  segment  centerline  and 
surfaces  are  defined  by  dx/dS  =  constant.  When  the  index  is  not  zero,  the 
segment  centerline  and  thickness  are  defined  with  fourth-degree  functions  of 
path  distance  from  the  transition  and  maximum  thickness  points,  respectively. 

The  specification  of  the  coefficients  for  these  functions  is  extra  input,  for  which 
the  format  is  shown  in  figure  12(c).  If  IDEF(IROW)  is  positive,  the  coefficients 
for  the  definition  polynomials  are  interpreted  to  be  functions  of  segment  length 
normalized  to  1.0;  but  if  lDEF(IROW)  is  negative,  the  coefficients  are 
interpreted  to  be  functions  of  segment  length  normalized  by  chord.  The 
reference  point  for  the  centerline  polynomials  can  be  either  the  transition  point 
or  the  segment  ends.  The  possible  combinations  are  shown  in  the  IDEF(IROW) 
summary  in  table  IV. 

designation  of  which  profile  loss  set  (I  variable  in  DLOS(K,J,l))  to  use  with  15 
particular  blade  row  .  If  the  input  value  of  lLOSS(IROW)  is  less  than  or  equal  to 
zero,  a  total  pressure  level  is  input  in  place  of  losses.  The  pressure  is  input  with 
the  parameters  shown  in  the  first  option  of  figure  12.  These  parameters  are 
stored  into  the  locations  of  loss  sets  4  and  5;  so  those  loss  sets  are  not  available 
for  use  with  any  blade  row. 

incidence  angle  (deg)  that  can  be  input  by  option.  If  the  tabular  option  is  used,  a  F10.4 


value  is  expected  for  each  streamline  starting  from  the  tip. 

molecular  weight  of  gas  (28.97  for  dry  air)  15 

number  of  annular  stations  at  which  radial  velocity  profiles  are  constructed  15 
during  computation 

number  of  blade  rows  (maximum  of  20)  15 

number  of  points  input  to  describe  hub  geometric  boundary  (maximum  of  40) 
number  of  loss  sets  input  (maximum  of  5)  15 

number  of  points  input  to  describe  tip  geometric  boundary  (maximum  of  40). 
number  of  sections  across  blade  for  which  fabrication  coordinates  are  desired.  If  110 


zero,  the  program  will  set  the  number  of  XCUT’s  on  the  basis  of  aspect  ratio. 

For  all  positive  values  the  program  will  set  appropriate  locations  to  represent  the 
blade.  Negative  values  of  NXCUT(IROW)  trigger  an  option  to  read  cards  for 
the  XCUT  values.  The  number  of  values  expected  for  a  blade  row  is  the  absolute 
value  of  NXCUT(IROW). 

number  of  streamlines  (maximum  of  1 1)  15 

option  controlling  amount  of  output  information  desired.  Interpretable  options  A4 
are  APPROX,  VEL.  DIA.,  DESIGN,  and  COORD.  If  the  first  four  characters 
input  in  OP  match  none  of  the  above,  the  program  will  try  to  proceed  with  the 
VEL.  DIA.  option.  The  program  completes  only  velocity  diagram  information 
when  run  with  the  APPROX  and  VEL.  DIA.  options.  With  the  APPROX 
option  the  locations  of  blade  edges  are  estimated  from  the  stacking  line,  but 
with  the  VEL.  DIA.  option  the  blade  edge  locations  are  input.  The  blade  edge 
data  are  read  from  extra  cards  at  the  end  of  the  data  set  for  a  particular  blade 
type.  The  axial  coordinates  are  temporarily  read  into  VTFI(I,J),  and  the  radial 
coordinates  are  temporarily  read  into  PO(l,J).  When  run  with  the  DESIGN  and 


Parameter 


Description 


Format 


OPM 


OPO 


PHI(IROW,J) 


PRA(IROW), 

PRB(IROW), 

PRC(IROW), 

PRD(IROW), 

PRE(IROW) 


PO(I,J) 


PR 


COORD  options,  the  program  designs  and  stacks  that  particular  blade  row. 

With  the  DESIGN  option  only  velocity  diagram  information  is  printed,  but  the 
blade  leading-  and  traiiing-edge  locations  are  for  the  stacked  blade.  The 
COORD  option  includes  the  printout  of  blade  section  properties  and 
coordinates  for  fabrication. 

additional  output  options  in  effect  if  OP  is  DESIGN  or  COORD  A4 


Card  column 

Additional  output 

7 

8 

9 

10 

O 

Off -design  punch 

T 

TSON1C  punch 

M 

Blade  elemeni  channel  microfilm 

M 

O 

M  and  O  options 

M 

T 

M  and  T  options 

_ 

_ 

_ 

.. 

_ 

Additional  output  options  in  effect  when  OP  is  COORD  A4 


Card  column 

Additional  output 

rn 

18 

20 

M 

Fabrication  coordinate  on  microfilm 

p 

Fabrication  coordinate  punch 

c 

MER1DL  punch 

M 

P 

M  and  P  options 

M 

C 

M  and  C  options 

ratio  of  inlet  segment  turning  to  outlet  segment  turning  (ratio  of  F10.4 

( dK/dS)]  \/\(dK/dS)2 1)  for  a  blade  element.  If  input  values  are  expected  by  use  of 
the  tabular  option,  the  data  cards  go  with  the  optional  cards  at  the  end  of  the 
data  set  for  each  blade  row.  A  value  is  expected  for  each  streamline  beginning 
from  the  tip. 

coefficients  for  polynomial  equation  to  define  profile  behind  b  row  P'  .jr.d  F10.4 
a  rotor  the  pressure  ratio  profile  is  specified  as 

p 

=  1.0  +  PRA*/?  +  PRB*/?2  +  PRC*/?3  +  PRD*/?4  +  PRE*/?5 

where  P,  is  the  stagnation  pressure  at  the  rotor  exit  tip  and  R  = 

(r,  - r)/(r,  - rh )— a  fraction  of  passage  height.  When  |PRA(IROW)j  >  100.0, 
another  option  is  activated.  The  input  profile  is  for  a  temperature  profile  T/T, 
instead  of  a  pressure  profile  P/P,.  The  data  value  of  PRA(IROW)  is  extracted 
from  the  input  value  by  adding  or  subtracting  100’s  until  the  remainder  is  in  the 
range  of  -  100.0  to  100.0.  At  a  stationary  blade  row  the  polynomial  is  for  the 
blade  row  exit  tangential  velocity  profile  in  ft/sec. 

V»  =  PRA//?2  +  PRB/R  +  PRC  +  PRD*/?  +  PRE*/?2  where  R  =  r/r, 

general  stagnation  pressure  array  in  Ib/ft 2  within  program.  The  I  index  is  the  F10.4 
station  index  and  J  is  the  streamline  index.  Only  (PO(l,J),  J  =  1 ,  NSTRM) 
values  are  input;  that  is,  the  streamline  value  for  the  first  calculation  station. 

The  input  values  are  read  in  units  of  psia. 

When  blade  edge  coordinates  are  input,  some  of  the  other  PO(l,J)  locations  are  F8.4 
used  for  temporary  storage  of  the  input  values  of  radius. 

desired  overall  compressor  pressure  ratio  FI 0.4 
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Parameter 


Description 


Format 


PTT(IROW), 

PTC(l.IROW), 

PTC(2,lROW), 

PTC(3,IROW), 

PTC(4,lROW), 

PTC(5,IROW) 

RHUB(l) 


ROT 

RTIP(I) 


SOLID(IROW) 

TALE(IROW), 

TBLE(IROW), 

TCLE(IROW), 

TDLE(IROW) 

TAMAX(IROW), 

TBMAX(IROW), 

TCMAX(IROW), 

TDMAX(IROW) 

TATE(IROW), 

TBTE(IROW), 

TCTE(IROW), 

TDTE(IROW) 

TILT(IROW) 


TITLE(l) 

TO(I.J) 


coefficients  that  describe  blade  row  exit  profile  when  it  is  input  as  an  option. 
PTT  is  the  blade  row  exit  pressure  in  psia  at  the  tip  (highest  radius).  The  other 
five  values  are  polynomial  coefficients  for 

P=PTT*(1.0  +  PTC1*R  +  PTC2‘R2  +PTC3-R3  +PTC4*/?4  +  PTC5*R5 

where  R  =  (r,  -r)/(r,  -rh)— fraction  of  passage  height  at  blade  row  exit 

radius  coordinates  of  a  set  of  points  that  define  geometric  hub  boundary 
(maximum  of  40) 

compressor  rotational  speed,  rpm 

radius  coordinates  of  set  of  points  that  define  geometric  tip  boundary  (maximum 
of  40) 

tip  solidity  of  a  blade  row  (ratio  of  chord  to  circumferential  spacing) 

polynomial  coefficients  of  ratio  of  blade  element  leading-edge  radius  to  chord, 
where  //*/c  =  TALE  +  TBLE*R  +TCLE*/?2  +TDLE*/?3 
where  R  (r,~r)/{r,  -rh)— fraction  of  passage  height  at  blade  leading  edge 


polynomial  coefficients  of  ratio  of  blade  element  maximum  thickness  to  chord, 
where  rmajf/c  =  TAMAX  +  TBMAX*R+TCMAX*R2  +TDMAX*R3 


polynomial  coefficients  of  ratio  of  blade  element  trailing-edge  radius  to  chord, 
where  /,e/c  =  TATE  +  TBTE*R  +  TCTE*/?2  +TCTE*/?3 
where  R  ( r ,  ~r)/(r,  —  r/, ) —  traction  of  passage  height  at  blade  trailing  edge 


angle  of  stacking  axis  tilt  (deg)  in  circumferential  direction  (r-8  plane).  The  angle 
is  positive  in  the  direction  of  rotor  rotation.  If  |TILT(lROW)i  >  100.0,  a  curved 
stacking  line  is  specified  according  to  r-rrej  =  C{ sin  7  -  sin  7 reA,  and  the  code 
of  the  TILT(IROW)  is— 


(x)xxxx-xx 


—  tilt  angle  at  hub  in  degrees  with  sign  of 
overall  TILT(IROW)  number 


—  tilt  angle  at  tip  in  degrees.  Circled  digit 
controls  sign  of  tip  tilt  angle.  Even  digit 
gives  tip  tilt  angle  same  sign  as  hub  tilt 
angle.  Odd  digit  gives  tip  tilt  angle 
opposite  sign  of  hub  tilt  angle. 


For  example:  12332.65  gives  a  hub  angle  of  23°  and  a  tip  angle  of  -32.65°. 

description  of  compressor  for  printout  and  later  identification 

general  stagnation  temperature  array  in  program.  Only  (TO(l,J),  J  =  l, 
NSTRM)  values  are  input;  that  is  the  streamline  value  for  the  first  calculation 
station.  The  input  values  are  in  units  of  °R. 


F10.4 

F10.4 

FlO.4 

F10.4 

FI0.4 

F10.4 

FI  0.4 

F10.4 


18A4 

F10.4 


Parameter  Description  Format 

TRANS(IROW.J)  location  of  transition  point  on  blade  element  centerline  as  fraction  of  blade  F10.4 
element  chord.  If  input  values  are  expected  by  use  of  the  tabular  option,  the 
data  cards  go  with  the  optional  cards  at  the  end  of  the  data  set  for  each  blade 
row.  A  value  is  expected  for  each  streamline  beginning  from  the  tip. 

VTH(I,J)  general  tangential  component  of  velocity  array  in  program.  Only  (VTH(1,J),  F10.4 

J  =  l,  NSTRM)  values  are  input;  that  is,  the  streamline  value  for  the  first 
calculation  station.  The  input  values  have  units  of  ft/sec. 

When  blade  edge  coordinates  are  input,  some  of  the  other  VTH(I, J)  locations  F8.4 
are  used  for  temporary  storage  of  the  axial  coordinates  of  the  points. 

XCUT(IC)  radial  location  of  blade  section  planes.  Whether  or  not  data  cards  are  read  for  F10.4 

values  of  XCUT(IC)  for  a  blade  row  is  controlled  by  the  value  of  NXCUT  (IC). 

Any  XCUT(IC)  cards  are  read  in  an  output  routine.  Therefore  they  must  follow 
all  cards  read  in  subroutine  INPUT;  that  is,  they  follow  the  ANNULAR  card 
for  the  last  calculation  station.  There  is  no  index  identifying  the  data  with  a 
particular  blade  row,  so  the  data  sets  for  the  blade  rows  are  expected  in  the 
order  that  one  would  see  the  blade  rows  in  moving  through  the  compressor  from 
the  inlet.  Start  the  set  of  points  for  each  blade  row  on  a  new  card.  It  is 
preferable,  but  not  necessary,  to  list  the  XCUT(iC)  for  a  blade  row  in  order 
starting  from  the  tip. 

XHUB(I)  axial  coordinates  of  set  of  points  that  define  geometric  hub  boundary.  The  axial  F10.4 

extent  of  the  coordinates  must  at  least  reach  the  first  and  last  calculation 
stations.  The  hub  coordinates  must  have  the  same  reference  origin  as  other 
input  axial  coordinates,  that  is,  casing,  blade  edge,  and  stacking  line 
coordinates.  The  number  of  points  input  should  be  4sns40. 

XTIP(I)  axial  coordinates  of  set  of  points  that  define  geometric  tip  boundary  (See  F10.4 

XHUB(I)  for  additional  comments.) 

ZHUB(I)  blade  data  set  hub-axial  coordinate.  When  the  data  set  is  a  blade  rather  than  an  F10.4 

ANNULAR  station,  ZHUB(I)  is  the  axial  location  of  the  blade  stacking  line  at 
the  hub. 

ZMAX(IROW.J)  location  of  maximum  thickness  point  as  fraction  of  blade  element  chord.  If  input  F10.4 
values  are  expected  by  use  of  the  tabular  options,  the  data  cards  go  with  the 
optional  cards  at  the  end  of  the  data  set  for  each  blade  row.  A  value  is  expected 
for  each  streamline  beginning  from  the  tip  with  a  leading-edge  or  transition- 
point  reference  according  to  option  (see  EB).  With  a  transition  point  reference 
the  values  input  are  (m  —  O/c 


ZTIP  blade  data  set  tip-axial  coordinate.  (See  ZHUB(I)  for  similar  additional  F10.4 

comments.) 
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Appendix  C 

Development  of  Equations  of  Motion  into  Form 
Used  in  Computer  Program 


In  the  computer  program  the  equations  of  motion 
are  applied  at  calculation  stations  that  are  presumed 
to  be  outside  the  blade  rows;  so  the  equations  of 
motion  are  more  conveniently  developed  in  an 
absolute,  rather  than  a  relative,  coordinate  system. 
The  general  equation  of  motion  (eq.  3(21)  of  ref.  7)  is 


in  terms  of  equations  (C3)  and  (C6)  the  term 
V  x  ( v  x  V)  can  be  expressed  as 


fx(7x  V)  = 


—  +  VH=  Vx(V  x  V)  +  tVs+f 

OT 


d(rVe) 

rdn 


When  steady  flow  is  assumed  and  the  local  friction 
force  is  ignored,  equation  (Cl)  reduces  to 

VH=Vx(Vx  V)  +  tVs  (C2) 

In  orthogonal  curvilinear  coordinates  the  velocity 
vector  can  be  expressed  as 

V=SVg+rhVm+nVn  (C3) 

where  m  is  in  the  streamline  direction  in  the 
meridional  plane  and  n  is  in  the  normal  direction  in 
the  meridional  plane.  Of  course  Vn  is  zero 
everywhere  for  this  application.  The  curl  term  in 
general  can  be  expressed  as 


+y  k  _  d^m  +  y  k 
dm  +  "  "  dn  +  mKn 


d(rVe) 

dVn 

,  n \dVm 

d(rVe) 

L  dn 

dd 

rl  dd 

dm 

where  km  and  k„  are  the  curvature  of  the  streamline 
and  the  normal,  respectively.  All  terms  containing 
Vn  are  zero  for  this  application.  The  assumption  of 
symmetric  flow  in  the  circumferential  direction 
makes  dVm/dd  equal  to  zero.  Also,  because  angular 
momentum  does  not  change  on  streamlines  outside 
the  blade  rows 


=  §!0]  +  m|0]+n  Y±*W} 
L  r  dn 


+  v  -  —  -V2  k 
T  m  dn  m  m 


Now  break  equation  (C2)  into  the  three  component 
equations.  In  the  6  direction 


dH  __  ds_ 
rdd  1  rdd 


The  zero  in  equation  (C8)  recognizes  circumferential 
symmetry  of  s.  In  the  meridional  plane  streamline 
direction 


dH_  jte 

dm  1  dm 


The  zero  in  equation  (C9)  comes  from  the 
assumption  that  entropy  does  not  change  along 
streamlines  that  are  outside  the  blade  rows.  In  the 
meridional  plane  normal  direction 


V»  d(rVe) 


+laS  (C10) 


*rv»)  -0 

dm 

Thus  equation  (C4)  reduces  to 


Equations  (C 8)  to  (CIO)  apply  to  the  three  curvilinear 
component  directions.  However,  in  the  program 
velocity  and  state  values  are  available  along  station 
lines;  so  it  is  of  computational  convenience  to  apply  a 
component  equation  along  a  station  line.  To 
accomplish  this  objective,  the  derivatives  in  the 
meridional  plane  are  converted  from  the  orthogonal 
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streamline  and  normal  directions  to  the  generally 
nonorthogonal  streamline  and  station  line  directions. 
The  angle  nomenclature  for  the  conversion  is  shown 
in  figure  13. 

The  enthalpy  gradient  in  the  station  line  direction 
can  be  expressed  as 


dH 

dt 


VHI 


Therefore 

dJ/Ve)  _  d(rV»)  1 
dn  dl  cos(a  +  X) 

dVm  dVm  /  ^  dVm  .  , 

dl  on  dm 


_  dH  d8  dH  dm  dH  dn 
~  rdd  dl  +  dm '  ~dl  +  Jn  '  di 

=  |0|-  [0]  +  |0|  sin(a  +  X)  +  ~  cos(a  +  X) 


Therefore 


dJjn=<Wm  1 
dn  d!  cos(a  +  X) 


dV, 


tan(a  +  X) 


dH 

dl 


— cos(o  +  X) 
on 


(Cll) 


ds  ds  _  , ,  ds  .  ,  . , 

-77  =  T-cos(a  +  X)  +  -r—  sin(a  +  X) 
dl  dn  dm 


In  general  a  station  line  derivative  can  be  expressed  as 


=  -t-  cos(a  +  X)  +  [0lsin(«  +  X) 
on 


d  d  dn  d  dm 
dl  dn  dl  +  dm  dl 


Therefore 


=  ~  cos(a  +  X)  +  ~sin(a  +  X) 
dn  dm 


(Cl  2) 


ds  ds  1 

dn  dl  cos(tt  +  X) 


(C13) 


(Cl  4) 


(C15) 


When  equation  (Cl 2)  is  applied  to  the  other  normal 
derivatives  of  equation  (CIO),  the  following  relation 
develops: 


d(rV6)  _  d(rVg) 
dl  dn 


cos(o  +  X)  + 


dirVj) 

dm 


sin(a  +  X) 


d(rVe) 

dn 


cos(a  +  X)  +  [0]sin(a  +  X) 


The  application  of  equations  (Cl2)  through  (C15) 
to  (CIO)  gives 


dH  _  Ve  d(rV$)  dVm 
dl  r  dl  m  dl 

~  y2mkmcos(a  +  \)  +  t 


■  dJj 

m  dm 

ds 

dl 


™  sin(a  +  X) 
(Cl  6) 


The  streamline  curvature  km  is 


k  —  l 
km~Jrn~R^ 


(C17) 


where  Rm  is  the  meridional  plane  streamline  radius 
of  curvature.  Substituting  equation  (CI7)  into  (C16) 
yields  the  following  form  for  the  meridional  velocity 
gradient: 


dVm 

dl 


dH  d(rV#) 

~di-v«-Tdr 


+  y m  -7^sin(a  +  X) 


V2  ,  xx  ds 
+  cos(a  +  X)  -  r  -T 

al 


(Cl  8) 
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Axial  direction 

figure  13.  -  Angle  nomenclature  for  direction  derivatives. 


The  state  properties  appearing  in  equation  (Cl 8)  are 
H ,  l,  and  s.  However,  two  state  properties  are 
sufficient  to  establish  the  others  at  a  point.  For  a 


thermally  perfect  gas  (p  =  p<Rt)  it  is  rather  easy  to 
compute  other  state  properties  from  two  selected 
properties;  so  it  is  desirable  from  a  computer  storage 
standpoint  to  store  only  two  properties  throughout 
the  flow  field.  The  two  properties  selected  were 
stagnation  temperature  and  pressure.  These  two 
properties,  along  with  the  velocity  components,  are 
sufficient  information  for  the  calculation  of  the  other 
state  properties.  If  these  two  properties  can  be  used 
directly  in  the  equations  of  motion,  the  need  to 
compute  some  state  properties  may  not  exist.  To 
express  s  in  terms  of  T and  P,  start  with  the  property 
relations 


—  =  dh  - t  ds 
p 


For  the  introduction  of  stagnation  properties  note 
that  the  thermodynamic  process  of  moving  between 
the  static  and  stagnation  states  is  isentropic  by 
definition.  Thus  equation  (C19)  for  this  process 
becomes 

&-dh 

P 

For  a  calorically  nonperfect  gas  this  becomes 


ds  _  1  dh  1 
dl~  t  dl  pi 


-  =cp(i)dt 
p 


-(fc> 


cp(t)dt 


dp  1  cp(t) 
p  'Git* 


[P  d-E  =-Lfr  cp{,)  dt 

Jp  p  (R  3/  t 


Equation  (Cl 9)  used  as  a  derivative  with  path 
distance  can  be  written  as 


ds  _  1  dh  _  1  dp 
dl~  t  dl  ~&di 


d^Pexp^-  1/(R|  cp(t)/t  dtjj 

X 

r  1  dh  1  dP  If  T  cD(t)  ] 

r  t  dl  Pt  dl  exp[  (R  3,  t  dtJ 

=  1^ -±dI(£.\+  P  (P\d\[TcP(t)A\ 

t  dl  pt  dt^P/  (Rpt '  p'  dt[)  i  t 

_l  dh  _<RdP  drfr  CjAt)  ] 

~ldi  Pdi+di[\,,dt\  (C22) 


The  application  of  Liebnitz’s  rule  to  the  last  term 
gives 


cp(T)  dT  cp(t )  dt 
T  dl  t  dl 


The  variable  (cp(t)/t)  is  not  a  direct  function  of  path 
distance;  it  is  a  function  of  temperature  alone. 
Therefore  the  partial  derivative  with  respect  to 
distance  must  be  zero.  Thus  the  derivative  of  the 
integral  can  be  expressed  in  terms  of  gradients  at  the 
limits  so  that 


r fr  Cp(Q  ^,1  cp{T)  dT  Cp(t)  dt 


l  <w_l  dh 
T  dl  t  dl 


Substituting  equation  (C20)  gives 


Substituting  (C23)  into  (C22)  gives 
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ds  _  1  dh  (R  dP  1  dH 
dl  t  dl  P  dl  +  T  dl 


ds=\_  M  _  ft  djP  =  l  dH  _  dj> 
dl  T  dl  P  dl  T  dl  p0  T  d! 


Equation  (C24)  is  essentially  equation  (C21) 
expressed  in  stagnation  state  variables.  Equation 
(C24)  would  turn  out  to  be  the  same  for  a  calorically 
perfect  gas.  Substituting  equation  (C24)  into  (Cl 8) 
gives 


„  dVm  dH  „  d(rVe)  ,  „  dVm ^ 

Vm~dT~~dT~  Ve~FdT  +  Vm  ~d^rsMa  +  X) 


,  .  . .  t  dH  ,  (Rt  dP 

1  cos(a  +  X)  -  —  ——  +  —  -JJ 


A  rearrangement  with  all  the  state  property  terms 
together  gives 


so  equation  (C26)  can  be  written  as 


1  Dh  „  n 
—  —  +  V  •  K=  0 
Dl 


Equation  (C27)  expanded  from  its  vector  form  is 


ifdh  VedJ,  dh_  dh\ 
a2\dt  r  M  +  ym  dm  +  n  dnj 


1  d(rVm)  1  dVe  1  d(rV„) 

r  dm  r  89  r  8n 


Outside  the  blade  rows  the  flow  is  assumed  to  be 
axisymmetric  and  steady.  Also,  because  there  is  no 
velocity  component  normal  to  the  streamline,  the 
equation  reduces  to 


'm  _(T-t\dH  l/a.dlnP  v<KrV9) 
7T  - v  -T" )  wr  +  ~  v<>  —rzr 


Vm  8h  s\  8(r  Vm) 
a2  dm  r  dm 


4"  V  m  km  —  0 


dV  V 2 

+  Vm  —in  sin(a  +  X)  +  cos(a  +  X) 
otn  a  f*\ 


Stagnation  enthalpy  is  defined  as 

y 2  i/2 

H  =  h+  -f-+Y 


All  the  terms  on  the  right  side  of  equation  (C25) 
can  be  computed  quite  accurately  except  dVm/dm, 
which  is  the  gradient  of  Vm  along  a  streamline  in  the 
meridional  plane.  The  distance  over  which  dVm/dm 
changes  sign  are  of  the  order  of  the  calculation 
station  spacing  so  that  representative  values  of 
dVm/dm  cannot  be  obtained  from  a  Vm  curve  fit 
along  meridional  streamlines.  A  better  value  of  this 
derivative  probably  can  be  obtained  by  means  of 
local  continuity.  From  equation  9(12)  of  reference  7 
differential  continuity  can  be  expressed  as 


dH=dh_.v  dVm  dV» 
dm  dm  m  dm  9  dm 


But  because  dH/dm  =  0  outside  the  blade  rows, 


dVm  V  dVf> 
dm  m  dm  9  dm 


Outside  the  blade  rows  angular  momentum  is 
conserved  along  streamlines;  so 


-  ^  +  V'M 
p  Dl 


d(rVe)_  dr  8Ve 

°-^T-J^y9  +  rd^ 


However, 


Rearrangement  gives 


_L  £!L 

q2  Dl 


dVe_  Ve  dr  _  Ve 

3 —  - -5—  = - sin  a 

dm  r  dm  r 
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Substituting  equation  (C31)  into  (C30)  gives 


Substituting  equation  (C34)  into  (C33)  gives 


Substituting  equation  (C32)  into  (C28)  gives 
a2  \  m  dm  r  /  r  dm 


m _  + 

-1  r 


+  -77  sec(a  +  X)  - 
dl 


tan(a  +  X) 


+  1^  +  Vmkm=0 


1“  "  OTTi  [H+Ot*  si""+  y"k")  <c3!> 

1  m 


The  curvature  of  the  streamline  normal  k„ ,  which  is 
da/drt,  needs  to  be  expressed  in  terms  that  can  be 
evaluated. 


da  da  ,  da  .  ,  .. 

^  =  ^cos(a  +  X)+^s.n(«  +  X) 


da  _da  1  da  si n(a  +  X) 

dn  dl  cos(a  +  X)  dm  cos(«  +  X) 


Calculation  of  dVm/dm  by  using  equation  (C35) 
should  give  a  somewhat  more  accurate  result  than  a 
curve  fit  or  a  finite  difference  computation  across 
increments  that  span  whole  blade  elements. 
However,  a  potential  divide-by-zero  complication 
has  been  introduced  with  the  term  M2m  - 1 .  In 
equation  (C35)  the  term  in  braces  in  essence 
represents  the  dA/A  term  of  one-dimensional  flow 
theory.  At  a  Mach  number  of  1.0,  dA/A  is  zero, 
which  is  the  throat  of  a  nozzle.  For  compressor  blade 
rows  the  throat  occurs  within  the  blade  passages. 
Internal  flows  adjust  around  locally  choked  regions 
so  that  the  throughflow  Mach  number  outside  the 
blade  only  approaches  1.  Computation  of  the 
detailed  nature  of  the  flow  is  not  available  from  only 
stations  outside  the  blade  row;  so  a  minimum  value  is 
imposed  on  the  denominator  through  an  empirical 
additive  term  to  help  stabilize  the  iterative  procedure. 
The  additive  center  term  is 


,  '<-1  ,n„2  , 

/=0.1  — — r exp  -  10 ,Mi,  -  1 

K-n 


,  da  da  tanfa  +  X) 

k"sdit  *  djsec(a  +  X)  fla¬ 


ils  characteristics  and  effect  on  the  denominator  are 
shown  in  table  V. 


Appendix  D 

Conic  Coordinates  of  Blade  Centerline  Path 


Local  blade  angle  is  defined  with  respect  to  the 
local  conic  ray  (fig.  14).  Let  the  blade  angle  vary  with 
path  distance  along  the  cone  according  to  the 
polynomial 

K  =  K,+aS  +  bS2+cSi  +dSA  (Dl) 

where  k,  is  the  blade  angle  at  the  transition  point 
between  segments  in  this  application.  The  path 
distance  S  is  with  respect  to  the  transition  point 
reference  but  always  positive  in  the  direction  from 
inlet  to  outlet. 

The  conic  radial  component  of  the  centerline  can 
be  found  by  integrating  the  differential  equation  for 
that  component 

dR  =  cos[*|dS  =  cos(xf  +aS  +  bS2  +  cS3  +  dS4  )dS 


Figure  14.  -  Blade  element  centerline  nomenclature. 


COS  K  =  1  - 


K2  K4  K®  K8 

2l  +  4T  "  6T  +  8T  '  '  ' 


(D3) 


When  equation  (Dl)  is  substituted,  the  terms  of  like 
powers  of  5  can  be  summed  to  give  in  symbolic  form 


(D2) 

The  problem  is  that  a  trigonometric  function  of  a 
polynomial  is  not  readily  integratable  in  closed  form. 
However,  the  function  can  be  expanded  in  series 
form  and  integrated  term  by  term.  Of  course  the 
series  is  infinite  but  it  is  convergent  within  the  range 
of  our  application.  In  the  following  presentation 
enough  development  is  given  to  show  the  form  of  the 
series.  Upon  application  in  the  program  a  tolerance  is 
used  so  that  no  more  terms  than  necessary  are 
calculated. 


cosk  =  [  ji +|  I2S  +  I  j3S2+!  4S3  -t- .  .  .  (D4) 


rs 


,  s2 

COSltds  =  (  ;IS+i  2  ~2 


S3 


S4 


+  1  ij  y  +1  k  -j  •  •  • 


When  terms  of  similar  coefficients  are  combined,  the 
following  form  evolves: 


cos  k  dS  = 


-  sin  k, 
a 


cos  K[  sin(aS)  +  sin  «,  cos(o5)J 

.  .  f  S2  a2  S5  a4  S7  06  S9  \ 

+  bS'n«t[-  y  +  TT"4!T  +  6!  9  +'"j 


+  b  cos  k,(^ 

S 4  a3  S 6 
-<7T  +  3!  T 

S8 
5!  8 

b2 

/  S3  a2  S1 

a4  59 

^  2  cos  k f 

F- 

1*^ 

1 

4!  T 

b2  . 

+  y  sin  K, 

(  S 6  a3  S8 

V  6  3!  8  ' 
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~  ***** 


,  .S1  a2d  S9  ad  S"  ad2  Sn  a S13 

+  b  cos  K(  ^  7  +  2  9  4!  11  +  2  2  +  6!  13 


s¥^  (  _ald  dl\?lL 

3! (2)  14  +\  8  +  3!  /  15  - 


.  .  a3d  S10  d 2  S“  ,  a3dS'2  a2d2  S13 

+  6 sin  Kt\ad  g  3,  J0  +  2  ,,  +  5!  12  2(2)  ,3 


a1  S14  a*d2  S15  /  <79cf  _  ad 3  \  S16 

7T  u  +  Sul)  iy+V'9r_ir/i6'. 


ft 2  .  rs9  a2d  Sn  a4d  S13  ad2  S14  a6d  S15 

~2  S,n  v  "9  2  FT  +  "4!  13  2  14  6!~  13" 

ft2  (  .S10  a3d  S12  d2  S13  a5d  S14  a2  d2 

"2  cos*t\ad-jQ  31  12  +  2  13  +  5!  14  2  (2) 

b3  /S11  a2d  S13  a4dSls  ad2  S'«  \ 

—  cosk,^  11  2  13  +  4!  15  2  16  +  "  7 

ft2  .  (  S'2  <J3d  S)4  $15  \ 

■  3!  S,n  K‘\-adl2  +  irT4-_Tl5-+"7 

b 4  .  /  ^S'3  fl2  S>3  \ 

4!  smK'V_dlT  +  TJT5'+--  7 


ft4 

+  -  cos  K, 


ft 5 

+  Jj-  COS  K, 


(-4-) 


+  csm  k, 


/  S4  ^ 

V  4  +  2  6  4!  8  +  '  ’ 


•) 


(  S5  a3  S1  a5  S9 
+  ccosk,^  a_5”"*’3!  7  5!  9 

c2  /  S1  a2  S9  \ 

+  ycosK,^-  y  +  y  -9  +  "  7 

c2  .  /  S»  \ 

+  y  sin  K,  y  + .  .  .J 

(  S8  a2tf  S10  a4d  S12  ad2  S13 
+ccos Kty—d-^  +  ~2  jq-  +  ^2"  +  ~2  HP  +  • 


+  csin  k, 


(  s9  a3d  S11  d2S>2  a5tfS>3  \ 

\act  9  3!  11  +  2  12  5!  13  7 


c2  •  a*S'3M  \ 

Ts,nKl(djr-Tdir+...J 
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+  .  .  . 


t 


C2  (  wS‘2  \ 

+  TcosK,^ir+...j 

.  .  (  Ss  a2  S1  a4  S9  \ 

+  dsin  + 

/  S6  a3  s8  \ 

+  rfcos«^-fl-  +  iyy+...j 

d2  (  S9  \ 

+  y  cos(C,(^-  -9  +•  •  -j 

fS“  a2  S13  ab  S14 

4-afccdcosK^  —  -  —  — 


a4  b2ac  S15 

ST  “  3!(4)  J  TS" 


a3  ft  be 

S'6 

_4!(2)  ~  T. 

16  +. 

a6  o3c  a2b2  c2 1 S17 
7f  +  4f(2)  +  (3!)2  “  IT  J  TT  ^ 


fc..  f  aS'2  bS"  ( a 3  r\S>4  a2fc  S13 

-after/ sin k, ^  2  12  2  ^  +^4,  2j  14  +3,(2)  ,5 


a5  aft2  a2c 
6!"  +  2(3!)  +  3l(2)J 


S16 

16 


a4b  abc  ft3 1  S17  ) 
5!(2)  IT  +  4!  JlT  j 


d2  .  (  S13  a2  S'1  \ 

+  abcYsin*l{-1T+JyTr+...) 

^  .  d2  (  a  S>‘  \ 

+abcJ  cos*,^  —  +...; 


With  these  groupings  shown,  patterns  of  terms  and 
coefficients  can  be  observed.  The  whole  equation  was 
coded  into  three  rather  brief  subroutines — one  for 
terms  with  two  coefficients,  COEF1  (two  of  the  four 
coefficients  a,  b,  c,  and  d);  another  for  terms  with 
three  coefficients,  COEF2;  and  one  for  terms  with  all 
four  coefficients,  COEF3.  Finally  the  coefficients  of 
the  terms  with  the  same  powers  of  S  are  summed;  so 
the  [  ]  terms  are  known  in 

*-*,  +  !  hS  +  lhy  +  [l,y 


•  +  [ 


S»_ 

n 


Because  in  the  following  developments  these 
coefficients  appear  frequently  within  parentheses, 
for  simplicity  the  [  )’s  are  replaced  with  c’s;  that  is. 


52  S3  S4 

R  —  R/  +C|S  +  C2  y  +  C3  y  +C4  —  +. 


Cn 


5! 

n 


(D6) 
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The  conic  angular  coordinate  can  be  expressed  as 


The  conversion  from  equation  (D6)  to  (D8)  begins 


e-t,= 


‘s  sin  k 
o  ~R~ 


dS 


(D7) 


where  both  sin  k  and  R  can  be  expressed  as  infinite, 
but  convergent  for  our  purposes,  polynomials  of  S. 
Since  a  polynomial  in  the  denominator  is  an 
undesirable  form  to  integrate,  the  polynomial  for  R 
was  converted  to  a  polynomial  in  the  numerator  of 
the  form  shown  in  equation  (D8). 


as 


R,  _  R, _ 

R  fl,+c,S  +  c2(S2/2)  +  cj(S3/3)+... 


_ 1 _ 

1  +  (C\/R,)S  +  (c2/R,)S2  +  (c3/R,)S3  +  .  .  . 


fs  sin  k 

’-■■Jo  irds 


![*«, 

R,  Jo  R 


sin  k  dS 


where 

^  =1  +RiS  +  R2S2  +R3S3  +  .. . 

A 


_  _ 1 _ 

1  -D,S-D2S2-DiSi-... 

where 


(D8) 


D  i 


£l 

3  R, 


,  etc. 


I  -£>,S-£>2S2  -O3S3  -• 


R 1 


R  2 


R  3 


r 


3  (  ^ 

I  +  0,S  +  (£>2  +  £>f)S2  +(C>3  +  2£»|£>2  +Oj)S3  +  <£>4  +  2£>,  £>3  +  D\  +  30,  £>2  +  D*t)S* 


1 

l-P|S-£>;S2  -P3S3 _ -D4S* 

0iS  +  02S2  +D3S3  +O4S4 

Di  S  -£)|S2  -D|D2Sj  -£>,D,S4 


(D2  +  £>f  )S2  +  (£)3  +  £>,  02)S3  +  (£>4  +  £>,  O3  IS4 

(P2  +  tf)S2  -  £>,(£>2  +  £>?)S3  -  P;(£>2  4  pfo4 

(£>3  +  2D,  £>2  +  D]  )S3  +  (P«  +D1D3  +D2  +  D2D\  )S4 
(£>3  +  2D,  D2  +  £>J)S3  -  P,  (£>3  +  2P3  +  2P,  D2  +  Pj)S4 

(P4  +  2D,  £>j  +  P^  +  3D?D2  +  P{  )S4 


✓  4 
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Table  VI  summarizes  the  preceding  division. 

The  coefficients  for  equation  (D8)  are  generated  in 
subroutine  RCOEF.  The  coding  for  the  procedure  is 
somewhat  complex,  but  in  general  not  much 
computation  is  required  to  satisfy  a  tolerance 
criterion  of  1.0E-08. 

The  conversion  of  sin  k,  where 

k  =  k,  -t  aS  +  bS 2  +  cS3  +  dS 4 

to  the  polynomial  form 

sin  k  =  /4 1  +-42S  +  >l3S2  +/I4S3  +/I5S4.  ■  •  (D9) 

is  accomplished  in  the  same  way  as  it  was  for  the 
cosine  series  (eqs.  (Dl)  to  (D5)).  In  fact,  the  cosine 
series  can  be  converted  to  the  sine  series  with  the 
following  substitutions: 


Cosine  series 

-sin  k, 
-cos  K, 
sin  k, 
cos  K, 


Sine  series 

cos  K 1 
—  sin  Kf 
cos  K, 
sin  k, 


1  [s  R,  . 

R  SmK 

1  fs 

=  jo(l+/?,S+/?2S2+/?3S2+.  .  .) 

x  (A  4  >4  2  5 4 A  2  S 2  4  /4  4  4  .  .  . ) 


*s;Io'4'+M2+iM',s 


4  (a4  3  4  /?  j  /4  2  +  /?2'4] 

4  (/4  4  +  R  \A$  4/?2^2  )S^  4  .  .  . 


1  l  .  „  .  A  2  +  R  i  A  i 


S2 


R'lA'S+  2 

,  A3  +RiAj  +  R2'4  1  C3 
+  -  S 


A4  +R\A 3  +R1A2  +  #3^1  ^4  + 


Consequently  the  same  routines  that  are  used  to 
compute  the  cosine  series  can  easily  be  modified  to 
compute  the  sine  series  coefficients  also. 

When  the  polynomial  series  coefficients  in 
equations  (D8)  and  (D9)  are  known,  the  integration 
for  £  is  straightforward. 


The  general  routine  for  establishing  the 
polynomial  coefficients  for  the  conic  coordinates  is 
EPSL2.  The  end  result  is  constant  polynomial 
coefficients  for  the  conic  coordinates  (R  and  e)  as  a 
function  of  S.  These  coefficients  are  saved  so  that  the 
conic  coordinate  at  any  S  of  interest  can  be  computed 
easily  with  subroutine  CONE. 
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TABLE  I.  -  OVERVIEW  OF  COMPUTER  PROGRAM 


Program  control 

Input  and  initialization 

Iteration 

Terminal  calculations 

Read  and  interpret  data 

Locate  calculation  stations 

At  each  station  for  each  streamline, 
estimate  stagnation  temperature 
and  pressure  and  axial  and  tan¬ 
gential  velocities 

Outer  loop; 

At  calculation  stations 

Set  coefficients  of  equation 

of  motion 

If  blade  design  option,  set 
incidence  and  deviation 

angles,  compute  new  blade 
edge  location,  and  reset 

calculation  station  location 

Inner  loop; 

At  each  calculation  station 

Solve  for  meridional  veloc¬ 
ity  distribution  to  satisfy 
equations  of  motion  and 
continuity 

Reset  streamline  location 

Overall  blade  row  performance 

on  streamlines  at  calculation 

station: 

General 

State  properties  (temper¬ 
ature  and  pressure) 

Velocity  diagrams 

Streamline  information 

Blade  rows 

Element  definition  param  - 

eters 

Incidence  and  deviation 

angles 

Aerodynamic  performance 
parameters 

Streamline  choke  margin 

Blade  section  parameters; 

Surface  coordinates 

Area,  moments,  etc. 
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TABLE  III. 


XIAL  LOCATION  HUB  AXIAL  LOCATION  TIP  BLOCKAGE  FACTOR  HUB  BLOCKAGE  FACTOR  MASS  BLEED  FRACTION 

(INCHES)  (INCHES) 


TIP  AXIAL  LOCATION  HUB  AXIAL  LOCATION  TIP  BLOCKAGE  FACTOR  HUB  BLOCKAGE  FACTOR  MASS  BLEED  FRACTION 

(INCHES)  (INCHES) 
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0000  9.3000  1. 0000  0.0000  0.5000 
0000  9  .  ft  0  0  0  1.0  0  0  0  0  0000  0.5000 
0000  11.1000  1.0000  0.0000  0.5000 
0000  16.0000  1.0000  0.0000  0.5000 
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63.26  6.79  0.00  45.93  45.86  40.02  33.88  9.69  1.1408  0.2946  0.7054  0.0288  0.2632  -0.248* 
81.31  7.53  -0.00  42.07  41.90  34.72  24.34  10.53  1.0300  0.2273  0.7727  0.0400  0.2715  -0.2S64 
98.07  7.88  0.00  36.90  36.96  28.28  11.35  11.29  0.8670  0-1531  0.7834  0.0773  0.3076  -0.3000 
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51  6!  ;  6>70  1.1417  3  6655  2.0252  7.0  91  1  1.5896  -  9.6  985  0.00  34  6.63  35.91  35.89  0.5504  0.2408 

61  70  1  6171  1  1692  3  6675  1.8648  6.529  9.8564  -9.4707  0.0038  7.52  27.81  27.85  0.5414  0.2934 
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3  3  1  6  37  0  1  1  5/  1  3  4  97  7  1.484  5  5.1  98  5.7  1  6  1  -8.5000  0.0046  1  0.39  1.51  2.31  0.5124  0.5058 

95  81  1.6164  1.1516  5  7616  1.2350  4.324  2.8865  -6.7263  0.0048  12.52  -19.10  -16.70  0.4856  0.6495 
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70.64  5.21  0.00  51.0*  51.05  42.64  40.36  12.19  1.1744  0.3831  0.6169  0.0522  0.0263  -0.0862 
82.47  5.47  0.00  50.39  50.33  39.86  35.05  14.10  1.1149  0.3405  0.6595  0,0765  0.0000  -0.1062 
97.12  5.12  -0.00  52.65  52.41  36.80  26.39  18.35  1.0280  0.2960  0.7032  0.1713  0.0000  -0.1615 
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97.36  1.6165  1.1577  2.0527  0.9204  5.566  5.3354  -7.6779  Q.Q148  12.40  -0.99  -0.19  0.4946  0.4105 
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